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ABSTRACT 



This thesis presents the results of an investigation of the lateral 
performance functions of an airplane by the analysis of the transient 
response of the airplane to finite pulse of control surface deflection. 

The method is an extension of pulse technique used by Clementson, 

Bef.(2), to measure the longitudinal performance functions of an airplane 
and further described by Seamans, Blasingmne, and Clementson, Ref. (3), for 
more general oases. In this thesis the successful extension to highly 
oscillatory lightly damped responses is shown. Also a possible technique 
of analysing divergent responses is indicated. 

The instrumentation and calibration neceseary to measure inputs of 
aileron and rudder deflection and outputs of time-rate of geometric yaw And 
roll are described in detail. The separate deflections of each aileron 
were added electrically to give the total differential aileron deflection 
thus obviating the necessity for adding these deflections during analysis. 
Microsyn signal generators were used to measure aileron and rudder deflections. 
Rate gyros were used to measure the angular velocities in roll and yaw. 

These were the only outputs measured due to the limited time available. 

Determination of the performance functions from the transient response 
was accomplished by approximating the Tourier Integral of the time responses 
of the inputs and outputs using the method of summing approximating triangles. 

The technique of measuring the transient response to finite pulse 
inputs is shown to be accurate, reliable and reproducible with a low level 
of uncertainties. 

Special methods of analyzing lightly damped oscillatory modes are 
developed and compared. The recommended method will vary somewhat with the 
form of the response curve but generally speaking the most satisfactory 
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method consists of approximating the Fourier Transform of the output over the 
early portion of the time response to a point beyond which the response Is, 
within the accuracy of measurement, only a damped sinusoid. At this value 
of time the Fourier Transform of the oscillatory mode obtained analytically 
is added vectorlally to the Fourier Transform of the initial portion of the 
response. 

Comparison with computed performance functions based on calculated 
derivatives shows qualitative agreement but clearly shows that these 
computed performance functions may differ considerably from the actual 
performance functions of the full scale airplane in flight. 

The saving in instrumentation and flight time effected by the pulse 
technique over that required by the sinusoidal technique is a major 
factor in recommending the former. Changes in the mass and moment of 
inertia due to fuel consumption are negligible for the short duration of 
test flights made possible by using the pulse technique. 
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CQBJJCT 



To Instrument a B— 25 J airplane and measure the 
lateral transient response* in rate of roll and rate of 
geometrio yarn to pulses of aileron and rudder deflection. 

To determine the lateral performance functions of 
the airplane from these responses by approximating the 
Courier Integrals of the Input an d output. 

To evaluate this method of obtaining these 
performance functions. 
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INTRODUCTIOH 



Recent demands for high-performance automatic control of aircraft 
for both Military and commercial usee hare focused attention on the 
airplane as the dynamic system being controlled. It is fundamental that 
the synthesis of a controlling system starts with the best possible 
evaluation of the dynamical system to be controlled. Additional interest 
in the dynamics of the airplane arises from the problems of stabilizing 
some of the auxiliary equipment utilized in modem airplanes, such as 
navigating and fire control equipment. 

The frequency response of an airplane can be calculated from computed 
derivatives based on wind tunnel tests. However, the comparison of these 
calculated responses with those measured on the actual airplane has shown 
appreciable discrepancies, due in part to the difference in stability 
derivatives obtained in the wind tunnel from those of the airplane, and in 
part to ignorance of the moments of inertia about the axes of motion. 

The most generally used methods of synthesis and design of automatic 
control systems requires the dynamic characteristics of the airplane to be 
presented as a frequency response. A measured frequency response is most 
often obtained by the direct method of applying sinusoidal inputs and 
measuring the amplitude and phase angle of the outputs. This is done for 
enough frequencies to define the complete frequency response. 

Pioneer work using the sinusoidal oscillation technique was carried 
out by the Cornell Aeronautical Laboratories with outstanding results . ^ 
These tests required an auto-pilot, a sinusoidal input generator and measur- 
ing equipment. Ixteneive flight time, with its consequent expense, is 
associated with this technique. 

An alternative technique, that of measuring the transient response to 
a step function input, has the drawback that such an input to the elevators 
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or aileron* will in a very *hort tine cause the airplane to rotate through 
a large angle and hence will hare displaced it to far from the original 
position of equilibrium that the motion i* non-linear. Such motion at the 
present time can he analyzed only by extremely tedious and lengthy method* 
of numerical computations . 

A technique which obriates a large part of the instrumentation and saves 
a great deal of flight time is that of measuring the transient response to 
a finite pulse of control surface deflection. This method was successfully 
used by Clement son^ Se ^ ^ in determining frequency response in the case of 
longitudinal motion. The mathematical background and the technique of 
application are presented in a paper by Seamans* Bias ingame, and 
Clement son . ^ 

In extending the pulse technique to the measurement of the lateral 

modes of motion of the airplane, the characteristic of spiral divergence, 

which exists in many airplane*, must be taken into account. The Fourier 

(Ref 3 ) 

Transformation method and it* approximation' * * are valid only for 
stable system*. Therefore, it becomes necessary to identify and remove 
the unstable mode prior to transformining the pulse response. Furthermore, 
the oscillation that occurs in lateral motion is lightly damped. This 
oscillatory mode must be given special treatment, or the analysis will 
become exceedingly tedious, in the regions of resonance. 

Methods of analyzing the oscillatory mode were successfully used in 
this investigation to find the lateral performance functions of the B— 25J 
airplane. However, this airplane did not have any measurable characteristic 
of spiral divergence f therefore, to illustrate a method of analysis when 
such a characteristic is involved, a simple system was analyzed. 
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or ailerons will in a very short tine cause the airplane to rotate throogh 
a large angle and hence will hare displaced it so far from the original 
position of equilibrium that the motion is non-linear. Such motion at the 
present time can be analyzed only by extremely tedious and lengthy methods 
of numerical computations. 

A technique which obviates a large part of the instrumentation and saves 

a great deal of flight time is that of measuring the transient response to 

a finite pulse of control surface deflection. This method was successfully 
(Ref 2 ) 

used by Clementson * * in determining frequency response in the case of 
longitudinal motion. The mathematical background and the technique of 
application are presented in a paper by Seamans, Bias ingame , and 
Clementson , ^ ) 

In extending the pulse technique to the measurement of the lateral 
modes of motion of the airplane, the characteristic of spiral divergence, 
which exists in many airplanes, must be taken into account. The Fourier 
Transformation method and its approximation^®^ are valid only for 
•table systems. Therefore, it becomes necessary to identify and remove 
the unstable mode prior to t ran s for mining the pulse response. Furthermore, 
the oscillation thAt occurs in lateral motion is lightly damped. This 
oscillatory mode must be given special treatment, or the analysis will 
become exceedingly tedious, in the regions of resonance. 

Methods of analyzing the oscillatory mode were successfully used in 
this investigation to find the lateral performance functions of the B-25J 
airplane. However, this airplane did not have any measurable characteristic 
of spiral divergence f therefore, to illustrate a method of analysis when 
such a characteristic is involved, a simple system was analyzed. 
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DXTIBMINATIOH OP AlEPLAHE LATERAL MOT I OH PI1P0BMAHCI TUNCTIOKS BT PUL SI 

TECHNIQUES 

The airplane in flight ma y he considered a* an operating component 
which operates on specified inputs to produce certain outputs. More than 
one output may result from a given input. The operation on an input to 
produce a specified output ie expressed mathematically by a concept of 
vide generality called a Performance Operator , 

By restricting the airplane and its controls to small excursions 
from an equilibrium condition, its motion may be defined by a system of 
linear differential equations. Then the response can be related to the 
input by means of a Performance Function . The concept of Performance 
Operator and Performance Junction were both developed by Or. C. S. Draper 
and his associates of the Instrumentation laboratory of the Massachusetts 
Institute of Technology. A complete discussion of their application to 
airplane dynamics is contained in Refs, (l) and (10) . The relation between 
the two may be summarized by the following equation: 



{PV - , + >‘ . , + > t (FT), , 



Vhen the output of an operating component can be related to the 
input by a set of linear differential equations, and the input is 
described by a complex exponential, the output is desoribed by 

j(v t + PA) 

^(out) * ^(ouOa # 



The ratio of the output to the input then defines the performance 



function of the operating component relating the two. In equation form 

J[v f t *7PA) out ] 



(pe) 






(ln)i 



J[ V + (pi) i»] 



The term e factors out, and it is seen that the performance 
function ie independent of time. Tor use in the synthesis of control 
or stabilising systems, the performance function is usually presented 
graphically, either as a locus on the complex plane or as an amplitude 
ratio and phase single plotted versus forcing angular frequency on the 
real plan*. (l)(2)(3)(l0) 



It is the need for the frequency spectrum of the airplsuie’t dynamics 
that has extended the study of these dynamics beyond the classical approach 
which, until the last few years, has been concerned with only the 
characteristics equation of the dynamic response and sensitivities . This 
approach does not give sufficient information to permit the design of a 
high-performance closed-loop system with the airplane as the controlled 
member* On the other hand, perform&noe functions provide this necessary 
information in a form readily adapted to the problems of closed-loop design. 

This investigation is concerned only with some of the performance 
functions which describe the lateral motion of the airplane. 

Specifically, they are: 

‘"Uv-V 

<PT) 4 V y 

(w) Jt» r .»j] 

The limitation in time available precluded the investigation of four 
more lateral motion performance functions of importance 9 e.g. 




«* 

T ] 



- 5 - 



(P/) A(. .» 

A 

(IT) A[. r .V ] 

(W) i£» r .« 

Analytical expressions for these performance functions can be obtained from 
the differential equations for lateral motion. 

The Laplace Transformations of the analytical equation* are: 

(* - Y y )v + Ur - p 0 

X XZ 

-tL + (- ^ s - L )r ♦ (• - L )p = L. 6 (*) 

T A XX r p a a 

I T7 

tH + (* - K )r ♦ (- s - H )p = H. fi (e) ♦ N. C (a) (l-l) 

T r l ll p 6 a a *r r 

Side force and rolling moment due to rudder deflection and side forcee 

due to rolling and yawing velocities are considered negligible. The axes 
and symbols conform to NACA notation. Fig. 1-1 picture* the axes, the 
symbol* are defined in the Glossary, Appendix I. 

In order to measure the outputs due to a specific input, all other 

inputs must be zero. Letting fi equal zero and specifying 6 a* the 

r a 

forcing function, the solution of equation (l-l) yields the performance 
function , 

'"'ivy 

where V is equal to r 

This performance function for the B-25J is computed in Appendix D, using 
values for the aerodynamic derivatives based on wind tunnel tests. Plots 
of the amplitude ratio and phase angle of this calculated performance 
function are given in Figs. 3-3 and 3-4, respectively. 

All of these lateral motion performance functions can be obtained 
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experimentally by sinusoidally oscillating the control surface and measur- 
ing the airplane’s steady state response. This sinusoidal method requires 
a sine generator and an auto-pilot to supply the inputs, and requires as 
well, the usual measuring equipment. In addition, considerable flight time 
is needed to obtain results over a sufficient range of forcing frequencies. 

In the case of rudder deflection inputs, a step-function input would 
yield a transient response from which the performance function could be 
determined. But a step-function of aileron deflection will 

produce a steady-state rolling velocity which in a short time will generate 
such a large angle of. bank that the motion could no longer be considered 
linear. 

The method of using a finite pulse of control deflection has the 
advantages associated with obtaining a transient response, and, since the 
control surface returns quickly to the trim position, the disadvantage of a 
step-function input is eliminated. Finite pulses of elevator deflection 
have been successfully used to determine the airplane’s longitudinal 
performance funct ions . the knowledge of the authors, 
this technique has not been applied to the lateral performance functions of 
an airplane . 

The mathematical background for a method of obtaining the performance 
function from the transient response to a pulse is presented in Refs. (2) and 
(3). Briefly, the method consists of approximating the Fourier Transforms 
of the input and the response. The ratio of the latter to the former is the 
performance function, that is, 




(TT)[q () (t)] 

(FT)[ qi (t)] 



( 1 - 2 ) 



or, more generally, 
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S = Jw 



(1-3) 



(FF) 






A[ Ql , q o ] ' (LT)[q ln (t)] 



Vfhen the response is convergent, it is shown in (2) and (3) that 

-JdAt w 

q (t )• 



(rr)[q o (t)] 



(FT)[UT(t)]. 0 2 



o n 



(1-4) 



n e 1,2,... 

The lateral response of an airplane requires special treatment "because 
of the lightly-damped oscillatory node and because of the divergent mode - 
found in many airplanes. The factored characteristic equation for the 
lateral motion of the B-25J, derived in Appendix D, shows both characteristics 

(1-5) 

The lateral response of an airplane having this characteristic equation 
would not be convergent. However, by resolving the response into components, 
the performance function can be determined i. 

For example, 

(PF), 



[s 2 + (2)(.14l)(l.34)e + (l.34) 2 ][e + 2.9][(s - .0013)] = 0 



( >KnVt> 



(LT) ^(out)(dlv.) (t) J 



(LT)[q lD U)] 



8 = JW 



(n)[q, 



(outHosc)^^ (rr)tq (out)(rem) (t)] 

* (FT)Cq^t)] * TmZjrn (1_6; 

The Laplace Transform of ^( ou t ) (div) ^ ^ *® °^^ a i ne ^ by first determining 

the analytical expression for ( 0Xi t ) (div) ^ • *^1® 1® obtained by plotting 

(on semilog coordinates) the response, q^ oQ ^)(t), for large values of time, 

where, for practical purposes, only the divergent mode remains. From this 

plot, the elope and intercept provide the necessary information to define 

the divergent mode analytically and the Laplace Transform of this analytical 

(ref 8 ) 

expression can be determined. * The Laplace Transform of the input is 

equal to the Fourier Transform, since the input has the value of zero for all 
values of time less than zero and its integral is convergent. Thus, the first 
term on the right side of equation (l-6) can be found from a time response. 

No evidence of a divergent mode was found in this investigation for the 
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B-25J . Therefore, it was unnecessary to employ the technique described in 
the above paragraph. However, to illustrate the procedure the performance 
function of an idealized physical system with divergence was determined from 
a computed time response. This performance function was compared with the 
ideal performance function, thus showing the feasibility of the method. This 
analysis is contained in Appendix C. 

The divergent component of the time response must be subtracted from the 
recorded response. The integral of the remaining function of time converges. 
However, because of the light ly—damped oscillatory mode, the approximation 
of the Fourier Transform by triangles becomes tediously long. This problem 
can be overcome by subtracting the oscillatory mode. The frequency, damping 
ratio, and amplitude of the oscillatory mode are determined as described in 
detail in Appendix D. Thi6 mode may be either subtracted over the region of 
time from zero to infinity, or it may be subtracted only in the region from 
any selected time, t^, to infinity. In the first case, the performance 
function becomes: 



w c (rT) Kout)(osc) (t) J * ^(out)(rem) (T 0^ (rT) ^ (UBT)(t) ^ 

( )[ V q J (FT)[UT(t)] £ 

1 n = 1,2,3 






+ EU (out)(rem) (T B >3. 

(rTHuTU)^ I U {ln )(T n )]<>' JnM i'' 



-Jn/W v 
o 



n ^ 1,2,... 

The unit triangle of t is here defined as a triangle of unit height and with 

a base length of 2At , as in Ref. 3. Similarly, the unit right triangle of 

t 1 s defined as a right triangle of unit height with a base length of At. 

Note that q, . w \(t)' is not the same function as q/ . w \(t) in eq. 
^(outMrem) ^(outHremJ 

(1-5). 

In the second of these two methods, the performance function becomes: 
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(FT)QUT(t )] S 

° n = 1,2, 




l • • • 



£ 



(1-7) 



(TT)[UT(t)] . I 

n = 1 




If the At associated with q, \ equals that associated with q , only the 



summations remain in the second term on the right side of equation (l-7). 



However, if Atq is not equal to At^, a correction factor must he applied, 
since the Fourier Transform of a unit triangle of t is 



The Fourier Transforms of both a unit triangle and a unit right triangle, as 
functions of wAt , are plotted in Ref. (3). The Fourier Transforms of both a 



developed in Appendix D, and their application to the determination of the 
lateral motion performance functions measured in this investigation is 
discussed in detail. 

Either of the above methods is considered acceptable. However, the 
results of this investigation indicate that adding of the oscillatory mode 
at tp is a somewhat better method than including the oscillatory mode from 
time equal to rero. 

The input pulse should have sufficient area to produce a desirable 
magnitude of response, but should not cause motion beyond the linear region 
of the airplane's aerodynamic characteristics. Also, to simplify the data 
reduction, it is desirable that the pulse be of short duration. The pulses 
used in this investigation were roughly estimated from previous data on 
static sensitivities, such as p b/2U, Ref. 4, together with the experience 



(sin wAt 

2 



(FT) (UT(t) ) * 
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gained from the longitudinal pulee teste on this airplane. The linearity 
ranges of the linkage* for the control deflection pick-off e (microeyn signal 
generators) and the sensitivities of the rate gyros were detemined on the 
basis of those estimates with a 100 percent safety factor. After flight 
trials, the rough estimate of the aileron pulse was increased by about 
50 percent In magnitude, with the duration remaining as originally 
estimated. On the other hand, the estimated rudder pulee was reduced 
slightly after flight trials. 
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FIG. 1-1. AIRPLANE REFERENCE AXES 



CHAPTER 2 



INSTRUMENTATION AND FLIGHT TEST PROCEDURE 
A. INSTRUMENTATION 

The airplane employed ae a test vehicle for this investigation was a 
US AT North American B-25J , Army Number 44-30328. The airplane had been 
stripped of all armament. A photograph, a three-view drawing and a table 
of general specifications and dimensions are presented in Appendix A. This 
airplane had previously been instrumented to record the normal acceleration 
and pitching velocity responses to a recorded pulse of elevator displacement. 

The amplifying and recording systems used in the earlier testing were retained 
for the present studies augmented by the following: 

1. Adjustable rudder stop device (mounted on the co-pilot’s right 
rudder pedal), fulfilling the dual purpose of providing rudder 
pulses of any desired fixed magnitude, and of locking the rudder 
while aileron pulses were being applied. See Figs. 4 and 5, 

Appendix B . 

2. Adjustable aileron stop device (mounted on the co-pilot's control 
column), fulfilling the dual purpose of providing aileron pulses 
of any desired fixed magnitude, and of locking the ailerons while 
rudder pulses were being applied. See Figs. 1, 2 and 3, Appendix B. 

3. Microsyn signal generator units, linked to each aileron for recording 
the magnitude of the aileron displacement. See Figs. 6 and 7, 

Appendix B . 

4. Microsyn signal generator unit mounted at the rudder control bell 
crank (right rudder) for recording the magnitude of the rudder 
displacement. See Fig. 8, Appendix B. 

5. Rate gyro, for measuring the angular velocity in yaw. See Fig. 15, 
Appendix B . 

6. Rate gyro, for measuring the angular velocity in roll. See Fig. 15, 
Appendix B . 

7. Electronic current regulator, used to regulate a set value of stiffness 
current applied to the elastic restraint generators of the roll and yaw 
gyros. 

Four channels of amplification and four channels of a Consolidated 
recording oscillograph were utilised to record: 
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1. Aileron input pulses. 

2. Rudder input pulses. 

3. Roll angular velocity. 

4. Taw angular velocity. 

The amplifiers and oscillograph are described in Appendix B. 

Several ground calibrations were made before any flights were 
attempted, and a final ground calibration was performed prior to the 
final flight. The methods of calibration and the calibration curves are 
presented in Appendix B. 

The center of gravity of the airplane was determined by means of a 
USA! electronic weighing kit. The center of gravity of the airplane was 
maintained near the center of gravity of the front main fuel tanks by 
shutting off the other tanks during flight. This caused the shift in the 
center of gravity during flight to be negligible. 

An accurate determination of the moments of inertia of the B-25J 
airplane for use in the theoretical calculations of this thesis posed a 
difficult problem. Obviously, they could be determined most accurately 
by experimental means, but, because of the size and weight of the B-25J, 
such a procedure was beyond the 6cope of the present study. As a practical 
alternative, the moments of inertia were fixed by careful estimates from 
known data Already established on similar aircraft. 

Estimates were first based on the fairly complete data for a B-25J 
contained in Ref. 1. However, a subsequent comparison of the values for 
lyy contained in the data with the value for 1^. determined for our test 
airplane by other investigators (Ref. 2) indicated that the test airplane 
and the reference airplane differed somewhat in configuration. Therefore, 
the final estimates were based on a comparison with similar data from the 
A—26 airplane discussed in Refs. 7 and 9. This comparison was expected to 
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increase the validity of the estimated figures because the test airplane 
and the A-26 had similar dimensions and differed only slightly in weight. 
This comparison yielded somewhat higher figures for 1^ and I z2 than those 
given in Ref. 1, in line with the higher value for I.j.-yd et ©rained experiment- 
ally for the test airplane, to which reference has already been made. 

B. FLIGHT TEST METHODS 

Thirty minutes prior to a flight, all equipment was turned On 
temporarily and the following items were checked to insure proper 
operation: 

1 . 300-volt d-c supply 

2. 26-volt three-phase gyro wheel excitation 

3. 150-volt d-c supply 

4. Gyro stiffness current 

5. Amplifier outputs 

6. Paper supply in oscillograph 

7. Frequency of the inverter output 

Soon after take-off, all equipment was turned on again this time 
for the duration of the flight — to insure sufficient warm up time, and 
the items mentioned above were rechecked. Then, with the airplane trimmed 
at 175 mph and 10,000 ft. altitude, test pulses were applied so that the 
zero position of the light traces of the oscillograph, the amplifier 
attenuation settings, and the input pulse amplitude could be adjusted. 

An effort was made to keep the attenuation settings constant so as to 
reduce the number of calibrations required. 

After the preliminary adjustments had been made, the airplane was 
carefully trimmed, with all controls held firmly against the stops, the 
pilot then started the oscillograph, and three seconds later gave the 
signal for the impulse. The oscillograph was stopped by the oscillograph 
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operator when the oscillations had attenuated to a negligible amplitude. 

The number of control stops used in these tests made it necessary for 
the pilot and co-pilot to follow a certain sequence of operations. This 
sequence was: 

1. Prior to recorded run: 

a. Pilot trimmed airplane to maintain straight and level flight, 
holding the elevator control against the elevator stop. 

b. Co-pilot adjusted aileron and rudder stops to give desired pulse 
amplitude*. (This was done by applying pulses of various 
magnitudes, within the linear range, until the desired response 
was observod on the oscillograph.) 

2. Tor a recorded run: 

a. Pilot trimmed airplane in straight and level flight, holding the 

, elevator control against the elevator stop. 

b. Co-pilot moved adjustable trim screw for rudder until the proper 
8 top was engaged f the rudder was then held firmly against this 
stop . 

c. Co-pilot moved the adjustable trim stop pin for the aileron until 
the proper stop was engaged, the aileron control was then held 
firmly against this stop. 

d. Pilot called oscillograph operator and, if the recording equip- 
ment was ready, pushed the button to start the oscillograph. 

a. Three seconds after starting the oscillograph, the pilot 
signalled the co-pilot to apply the pulse. 

On each flight, two records were made for each of the following types 
of pulse: 



1. Positive aileron (right aileron down) 

2. Negative aileron 

3. Right rudder 

4. Left rudder 

After the abotre records had been taken, the airplane was held on a 
straight course at 10,000 ft. for a sufficient time to make calibration 
records. The method of calibration is discussed in detail in Appendix B. 
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In addition to the oscillograph records, the following data were also 



recorded for each run: 

1. Attenuation settings 

2. Direction of control deflection 

3. Altitude 

4. Indicated airspeed 

5. Outside air tecperature 

6 . Fuel load 

7 . Engine rpm 

8. Inverter frequency 

9. Gyro elastic restraint current 

10. Signal generator excitation current 
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CHAPTER 3 



RESULTS AND DISCUSSION 



A satisfactory instrumentation and calibration of the test vehicle 
was achieved. It is shown in Appendix B that the maximum uncertainty in 
calibration was '1.7 percent. 

The flight test data was reduced by the methods outlined in Ref. 3 
with the following results: 

(a) The uncertainty in reproducing results, obtained by analyzing 
three independent records of rolling velocity response to input aileron 
pulees, was found to be 1.5 percent average deviation from the average 
measured performance function in amplitude ratio, and 2° average deviation 
in phase angle. (See Figures 3-1 and 3-2). 

(b) The comparison of the experimentally determined performance 
function with the performance function calculated from wind tunnel stability 
derivatives for rolling velocity response to aileron input showed very close 
agreement in phase angle at all frequencies within the range investigated 

(0 .5 to 6.0 rad/sec). The comparison in AR showed both curves following 
the same general pattern with the calculated performance function having 
a lower AR at all frequencies except the resonant frequency. The resonant 
frequencies were displaced both in amplitude ratio and in frequency in 
accordance with the respective differences in damping ratio and undamped 
natural frequency of the oscillatory modes. (See Figures 3-4 and 3-2). 

(c) The sensitivities of the calculated and experimental performance 
functions for rolling velocity response to aileron input compared as 
follows: 



S (A)(« O(calc) l/ ”° 

a a 

S (A)U V), w --1-00 

a X (exper ) (aver ) 
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Utilizing vb 8 data obtained from flight test data (Ref. 4) an 

cU & 

additional sensitivity is available, i.e., 



S ( a) ( fi V ) = " 0,9 ° 1 / eec 

a X tflt. test data) 

(This sensitivity allows a rough check on overall calibration of the 
instrumentation system) . 

(d) Comparison of the experimental performance function for yawing 
rate out to aileron input with the calculated performance function shows 
good correlation in AR. The PA comparison is satisfactory for low 
frequencies (.5 to 3 rade/sec). However, at high frequency a wide scattering 
of points is produced using three different systems of analysis. (See 
Figures 3-3 and S-4). 

(e) Sensitivity comparison of the calculated and experimental perform- 
ance functions for yawing velocity response to aileron input is as follows: 

Vcalc ‘ ' 209 1/>eC 

S < A >E‘ a "zbxp.r * l/,eC 

(f) From an investigation of three different methods of dividing the 

output into comnonents, as discussed in Appendix D, it was found that the 

most accurate and fastest system of graphical analysis for use in the case 

of a lightly damped oscillatory response was method (2) of Table 2, 

Appendix D. In this method the usual triangular approximation is used 

until the time response is only a pure damped sinusoid. At this point, 

t = t , the Fourier Transform of a cosine or sine, whichever is appropriate, 

P 

is added. 

(g) Investigation of the scattering of points at high frequency in 

the performance function, (PF)^^ y )• showed that the uncertainty was 

a Z 

a function of the following: 
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(1) The measurement uncertainty in determining rector magnitudes from 
the oscillograph record. 

(2) Determination of the proper base line, 

(3) The time response pattern — depends on whether the maximum ordinate 
occurs early or late in the tirte response, and on how large a proportion of 
the responee is purely oscillatory. 

(4) The ratio of magnitudes of the last plotted rector in the output 
summation to the resultant vector representing the response. 

(5) The method of component breakdown. 

(6) The frequency — uncertainty increases with frequency in the usual 

case . 

(7) Site of At used in the approximation of the input and output. 

Instrumentation and the actual test flights for this project consumed a 
relatively large percentage of the time available. Consequently, the time re- 
maining for analysis of the oscillograph records required that only those 
analyses most pertinent to the objectives of the project be undertaken. Since 
the main objective of the project was to investigate the practical aspects of 
the application of the pulse technique to dynamic systems in general, and since 
determination of the performance functions of the test vehicle were of secondary 
importance, the scope of data analysis was limited to the following investigation 

1. Determination of the uncertainty in repeating results by obtaining 

a performance function from each of three oscillograph records independently, 
using aileron pulse inputs and rolling velocity response. 

2. Ivaluation of several different graphical plotting methods as to 
accuracy, speed, and applicability to various response configurations. 

3. Comperisoh of the experimental performance function with the 
calculated performance functions, using response in rolling velocity and 
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yawing velocity to aileron inputs as examples. 

4. Seeking out of limitations imposed on the accuracy of the pulse 
technique due to the approximations involved in the process of converting 
from the time domain to the frequency domain. 

Table 1 of Appendix D lists the oscillograph records analysed, and the 
method of analysis used. 

The uncertainty in reproducing results using the pulse technique is 
considered to he entirely satisfactory, and well within the demands of usual 
engineering requirements. In evaluation of the ability of the pulse 
technique to reproduce results two investigations were pertinent. They 
were: (l) the ebility to produce the same performance function by analyzing 

different records, and (2) the ability of two different analysis methods to 
produce the same performance function when applied to the same oscillograph 
record. The results of the first investigation are shown in Figures 3-1 and 
3-2. The three oscillograph records used to obtain the experimental results 
shown in Figures 8-1 and 3-2 are tabulated in Table 1, Appendix D. The 
experimental curve shown is the average of three performance functions deter- 
mined independently. The maximum AH deviation from the mean curve is 5.9 
percent at u>^, = 1.5 rad/sec., and the average AR deviation from the mean 
curve i 8 1.5 percent. For the phase angle, the maximum deviation is 6° 
and the average deviation is 2°. 

The uncertainty in repeating results by utilizing the three graphical 
methods listed in Table 2 of Appendix D, was found to be a variable depending 
on the response pattern in the time domain, and the forcing frequency. In 
general, the method described as Method (2) in Table 2, Appendix D, offers 
the least uncertainty. This method utilizes the same degree of approximation 
in determining the Fourier Transform of the input and output from the records 
in the time domain out to t = t . Disregarding frequency effects, the 
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uncertainty of this method may be held fairly constant with variation in 



response pattern by selecting t at increasingly higher values as the 

P 

response pattern becomes increasingly more oscillatory. Comparing the 

pattern of rolling velocity response with the pattern of yawing velocity 

response using aileron inputs for both, it was found that in the former 

a low proportion of the response was oscillatory while in the latter a high 

proportion of the response was oscillatory. Consequently, in the latter 

case, if t was not increased beyond the t used in the rolling response 
P P 

case, a higher proportion of the response was represented by a Fourier 
Transform which is exact, while the input remains an approximation to the 
exact Fourier Transform. The ultimate limit is to let t Increase to a 



P 

time beyond which the response is negligible, plotting all vectors directly 
from the record. This is lieted in Table 2, Appendix D as method (l) . 

This method offers the least error in the hypothetical case where 
instrumentation inaccuracies and noise are at such a low level that 
ordinates may be measured accurately at large values of time. 

Increasing t as the percentage of oscillatory portion of the response 
increases is also important in nullifying errors occassioned by choosing t 

P 

at a point which i6 not exactly a true peak or zero of the response curve. 

An error of this nature serves to shift the phase angle of the added vector 

for the oscillatory remainder (t « t ) by an amount equal to u> .At . If the 

P ' 

magnitude of the last vector plotted in the summation of triangles is only 

a small percentage of the added vector, the error in phase angle may be 

large. However, if the magnitude of the vector added is of exactly the same 

length as the last vector of the summation, an error in t as large as t » hr 

P 

will yield no error in the PA of the response vector. 

The method of dividing the output into components, described in Appendix 
D, Table 2, as method (3), appears to be subject to the greatest inaccuracies 
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any of the methods Investigated in analyzing a lightly damped oscillatory 
response. Theoretically, this method of analysis is as sound as any of the 
methods investigated. However, the practical application of this method 
presents an additional source of error since the input and output are not 
represented with the same degree of approximation. In the lightly damped, 
highly oscillatory response, essentially the entire output is represented 
by the Fourier Transform of the oscillatory mode, which is the exact 
transform of an analytical expression. On the other hand, the input is 
approximated by a vector addition, the exactness of which depends on the 
accuracy of curve fitting achieved. At relatively high frequency, the 
small resultant vector of the remainder summation, although approximated 
with the same degree of exactness as the input resultant vector, is 
ineffectual in changing the total output phase angle which is dominated by 
the large vector representing the oscillatory mode. Consequently, the output 
phase angle becomes nearly static as the limiting phase angle of the oscillatory 
mode is approached, while the input phase angle continues to grow as the 
forcing frequency increases. 

In performing the component breakdown necessary in this system, additional 
measurement errors may accumulate if more than one set of measurements are 
taken from the response. That is, in removing the oscillatory mode, the 
ordinates as determined by the analytic expression of the oscillatory mode 
should be subtracted algebraically rather than graphically from measured 
ordinates of the total response. Furthermore, since the oscillatory mode 
represents the major part of the response, great care must be exercised in 
determining W^, DR, and 

In Figures 3-3 and 3-4, the results of analyses by the various methods 
discussed are shown. The scatter of points at high frequency shows the 
inaccuracies of the several methods. Theoretically, the PA curve should 
level off in the vicinity of -90° phase angle. However, Fig. 3-4 shows the 
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experimentally determined PA lagging lees at of 5 and 6 rad/sec than 
at 3 rad/eec. In investigating this discrepancy between theoretical and 
experimental results, the question of accuracy in approximating the Fourier 
Transform arose. 

To determine how closely the approximation of input and output time 
response by summation of triangles fits the actual curres, frequency spectra 
of yaw response and aileron input were prepared by numerical means. The 
results are shown in Figures 3^5 and 3-6. The circled points show the 
results of the numerical analysis and the dotted. curre shows the mean 
curre of the experimental points obtained in Figures 3-3 and 3-4. The phase 
angles determined by the numerical analysis show a variation of more than 30° 
from the average curve at many points. The explanation for such a discrepancy 
was found to lie in a difference in measurement of the ordinates, [q(t Q )] , 
between those used in the numerical analysis and those used in the graphical 
analysis. The measurements of the respective areas, for each analysis, are 
as follow: 







ARIA 


APTCA 






Numerical 


Graphical 






Analysis 


Analysis 


fit 


positive lobe 


0 .268 in sec 


0.26 in sec 


8 1 


negative lobe 


0.6170 in eec 


0.5990 in sec 


ad 


positive lobe 


0.1126 in sec 


0.1096 in sec 



In the two most important lobes, the deviation in measurement of area 
is approximately 3 percent, and in each case the graphical analysis is 
ths smaller area. Thus it is noted in Figure 3-7 that a small uncertainty 

of 3 percent in area determination may cause ae much as 30° difference in 

/ 

phase angle, provided . the vector plot of the output ’ ' winds up * ' near the 
origin. Since the numerical analysis gives only a check on how closely 
the triangular summation approximates the area of the curve, the performance 
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function was again determined at several forcing frequencies using the same 
ordinates as used in the numerical analysis. The triangles in Figure 3-6 
show the points thus determined. It is seen that the uncertainty of repro- 
ducing results is satisfactory, indicating that the use of triangles in 
approximating the area yields approximately the same accuracy as the numerical 
approximation. It is now apparent that any inaccuracies in measuring the 
ordinates hire % critical effect on the performance function when the end 
point of the vector summation is near the origin. However, the PA variations 
due to inaccuracies in measurement of ordinates is not constant, "but is a 
function of frequency. Although the measuring error in this case gives an 
increase of ordinate magnitude in both the positive and negative lobes, the 
errors in a positive lobe do not nullify the errors in a negative lobe 
except at certain frequencies which are indeterminate unless the actual 
measurement errors are known. Figure 3-7 shows the results of measurement 
error at u>^ ■ 5 rad/sec. Due to the phase angle shift with each added 
vector (u»^,At) the errors are apportioned in accordAnce with the sine and 
cosine of the vector phase angle multiplied by the magnitude of the error 
in that vector. If the sum of the imaginaries and reals of the measurement 
errors add up to give a vector which has an angle corresponding to that of the 
original resultant vector, no change in PA will result. However, this rarely 
happens, and at high frequencies a large phase angle change results. Hote in 

Fig. 3-7 that the largest differences in the vectors of the first two lobes 

o 

occur in measuring vectors which have phase angles close to +90 , i.e., vector 
10 in the first lobe and vectors 28 and 30 in the second lobe. All of these 
difference tend to increase the phase angle at • B. Eaoh frequency will 
have a different total PA error due to measurement uncertainties, since the 
phase angle of each vector changes with each different forcing frequency. 
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Again note that measurement errors hare an insignif leant effect unleae 
the end point of the summation ie near the origin. At e 1.75 rad/eec 
the PA Is the eame, whether determined graphically or numerically, and the 
effect of measurement differences on uncertainty in determining the 
performance function is negligible as shown by the superimposed points 
lying on the average experimental curve. 

In general, it appears that any uncertainty, whether in measurement 
of ordinates or base line location, gives a magnified uncertainty in PA 
if the end point of the summation is near the origin. 

The phase angle uncertainty ie not due to a breakdown in the theory of 
the pulse technique, nor caused directly by the triangular approximation 
of the input and response. Bather it ie a function of uncertainty in 
measurement of ordinates and in plotting which accumulate under the 
conditions of high frequency (rapid phase angle change of each added 
vector) and low amplitude ratio (requires and point to be near the origin), 
furthermore, this uncertainty in phase angle is not restricted to the 
predominately oscillatory response. It merely appears at a lower frequency 
in this case due to the rapid attenuation in amplitude ratio (small 
vectors in the early portion of the response) . The same phase angle 
uncertainty would occur in the rolling velocity reponse to aileron 
input if the performance function determination had been carried out to 
those high frequencies where the amplitude ratio was greatly attenuated. 



(A R>a a .w„,ir> 5W|£c 




^IG. 3-i. AR OF AIRCRAFT PERFORMANCE FUNCTION 
FOR ROLLING VELOCITY RESPONSE TO AILERON INPUT 

[PF](A)(£ Q ,w x ) 
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A 




FORCING FREQUENCY, Wf, in RAD/SEC 



FIG. 3-Z. PHASE ANGLE OF AIRCRAFT PERFORMANCE 
FUNCTION FOR ROLLING VELOCITY RESPONSE TO 

AILERON INPUT, (PF\^ q)W)c )- 
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AMPLITUDE RATIO, AR, in DEO PER SEC/DEG 




METHOD of ANALYSIS (App. D -Table 2} 
Method Symbol 

2. + Cosins added at tp 

l A Sine added at tp 

3 0 Rt. Tri + Osc. + Rem. 



FIG. 3-3. AR OF AIRCRAFT PERFORMANCE FUNCTION 
FOR YAWING VELOCITY RESPONSE TO AILERON INPUT, 

(PF’WAi .w*) 
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PHASE ANGLE, PA, in DEGREES 




FORCING FREQUENCY, Wf,i.n RAO/SEC 



METHOD of ANALYSIS (App. D -Tabl# 2) 
Method Symbol 

2 j + Cosine added at t p 

' a Sine added at t p 

3 O Rt.Tri. * Osc. * Rem. 



FIG. 3-4. PA OF AIRCRAFT PERFORMANCE FUNCTION FOR 
YAWING VELOCITY RESPONSE TO AILERON INPUT, 
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AMPLITUDE RATIO, AR, in DEC ppr SEC/ DEG 



O FREQUENCY SPECTRUM BY 
NUMERICAL ANALYSIS +SINE 
WAVE ADOED AT tp 




FORCING FREQUENCY, W{, in RAD/SEC 



FIG. 3-5. AR OF AIRCRAFT PERFORMANCE FUNCTION, 
YAWING VELOCITY RESPONSE TO AILERON INPUT, 

A)(cS a ,W z ) 
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CHAPTER 4 



CONCLUSIONS AND RECOMMENDATIONS 



A. CONCLUSIONS 



As a result of this investigation the following conclusions are 
submitted : 

1. The pulse technique is a practical and suitably accurate method 
for determining aircraft performance functions from flight teets. 

2. Economy in flight time and instrumentation are gained by the 

i 

pulse technique as compared to the einusoidal method. 

3. The accuracy of the re suite obtained is enhanced by the short 
duration of test flights in that the changes in mass and moments of 
inertia due to fuel consumption are negligible. 

4. Erroneous performance function phase angles will in general be 
obtained at relatively high frequencies when the amplitude ratio has 
attenuated to a very small value (order of five percent). This is due to 
the fact that the errors in measurements of the ordinates of the recorded 
response accumulate at those frequencies where the series of vectors, 
representing the ordinates and their individual phase angles, encircle 
the origin. 

5. Lightly damped oscillatory responses may be analyzed with reasonable 

simplicity and accuracy. However, in the case of a response which is almost 
wholly oscillatory, uncertainty in the data may cause anomalous phase effects 
to appear in the plots of the performance function at the higher forcing 
frequencies when the amplitude ratio is considerably attenuated. These phase 
effects are more pronounced in the performance functions which are highly 
dependent on coupling terms in the equations of motion rather than the 
forcing function itself, since high attenuation in AR occurs at a relatively 
low frequency. For example, PTr^ w and PFj- fi are such performance 

- L a ZJ L r » X J 

funct ions • 
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6. A predominately oscillatory response is best approximated by use 

of isosceles triangles to a value of time beyond which only a damped sinusoid 
remains. The Tourier Transform of the sinusoid is then obtained by trans- 
forming the analytic expression which is then added vectorially to the 
approximate transform of the early portion of the response. 

7. The most important sources of error in reduction of data are: 

a. Choosing the base line on the oscillograph record. 

b. Measurement of the ordinates from this base line. 

8. The calculated and experimental performance functions are 
qualitatively comparable, but the quantitative differences indicate the 
requirement that the dynamic characteristics of all high speed and non- 
conventlonal aircraft be determined by flight tests. 

B . RECOMMENDATIONS 

To aid further research in the use of the pulse technique to determine 
aircraft performance functions the following recommendations are presented: 

1. Obtain relatively long records of lightly damped oscillatory 

responses to facilitate the determination of W and DR of the oscillation. 

n 

2. Amplify low amplitude predominately oscillatory responees as much 
as is consistent with the noise level and linearity of the recording 
instruments, the friction torque level of the gyros, and the linearity of 
aircraft motion. 

3. Limit the shock mounting of the gyros to landing and take-off 
only. During tests the gyros should be rigidly attached to the aircraft 
structure so that the input to the gyros is the true angular velocity of 
the Airplane and is not modified by the flexibility of the shock mounts. 

4. Superimpose a reference trace on each response to aid in the 
determination of the base line and in indicating any drift or divergence 
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of the response trace. 

5. Investigate by flight teet the lateral notion performance 
function# of an airplane haring a relatively high degree of spiral 
divergence. Since the degree of spiral divergence increases with angle 
of attack, it nay be necessary to perfom these measurements at a high 
angle of attack in order to obtain suitable responsee. 

6. Investigate the performance functions of an airplane having 
dynamics quite different from those of the B—25J . Several modem fighter 
types exhibit an undamped oscillation about the Z axle. A study of this 
motion should prove most worthwhile and interesting. 
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APPENDIX A 



DESCRIPTION OF AIRPLANE IN TEST CONFIGURATION 
The airplane used in the tests was a USAF B-25J under bailment contract 
to the Instrumentation Laboratory, M.I.T. The airplane had previously been 
stripped of all guns and turrets, as shown in Fig. A-l. 

Following is a listing of the specifications and dimensions of the air- 
plane as it was used in the tests described in this thesis: 



Manufacturer 


North American Aviation 
Corporation, Inc. 


Type 


R-25J-25-N.C. 


Serial Number 
Center-of-gravity range 


44-30328 


Forward limit 


16.5 percent M.A.C. 


Aft limit 


34.0 percent M.A.C. 


For tests 


27.4 percent M.A.C. 


Gear retraction moment/1000 


59 .7 in- lbs 


Overall length 


53 ft 5.75 in 


He ight 

Wing, airfoil section 


16 ft 4.19 in 


Root 


NACA 23017 


Tip 


NACA 4409 R 


Span 


67 ft 6.704 in 


Area (Total) 


609.8 §q ft 


Area (less ailerons) 


577.67 sq ft 


Aspect ratio 


7.48 


Taper ratio (tip chord/root chord) 


0.415 


Chord at root 


12 ft 10.6 in 



Chord near tip (380.75 from fuselage) 6 ft 10.14 in 
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Projected tip chord 
Mom aerodynamic chord 



5 ft 4.257 in 



Length 

Dietance of leading edge "back of 


116.16 In 


noee reference datum line 


212.68 in 


Incidence 

Hoot 


3° 0* 30’ ’ 


Tip 

Dihedral (25 percent line) 


0° 21’ 39” 


Wing center section 


4° 38’ 23” 


Wing outer section 


0° 12’ 39” 


Sveepback (leading edge) 
Aileron 


4° 12’ 13” 


Type 


Internal balance unsealed 


Area 


32.13 sq ft 


Span 


137.25 in 


Chord 

Inboard 

Outboard 


19.87 in 
12.45 in 


Travel 

Up (from neutral) 
Down (from neutral) 


o c 

3S 


Tab area 


3.86 sq f t 


Tab span 


57.911 in 


Tab chord (mean) 


5.375 in 


Tab travel 

Up (from aileron trailing edge) 
Down (from aileron trailing edge) 

Wing flaps 


c o 

22 


Type 


Yariable-slot 
trailing edge 


Area (total) 


75.8 sq ft 
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Center section flap* 



Area 


32,6 sq ft 


Span (mean) 


127.5 in 


Travel (down) 
Outboard flap 


43° 


Area 


43.2 sq ft 


Span (mean) 


214.5 in 


Travel (down) 
Horizontal tail surface 


o 

to 


Area (including elevator*) 


132.4 sq ft 


Span 


22 ft 2 in 


Maximum chord 


7 ft 1.875 in 


Incidence 


2° 


Dihedral 


None 


Ilevator area (two, including tabs) 


50.6 eq ft 


Ilevator spaa (two) 


216.5 in 


Ilevator travel 

Up (from streamline 
with stabilizer) 


26° 


Down (from streamline 
with stabilizer) 


10° 


Ilevator trim tab area (total) 


4.22 eq ft 


Ilevator trim tab span (2) 


57.2 in 


Ilevator trim tab mean chord 


7.875 in 


Ilevator trim tab travel 

Up (from elevator trailing edge) 
Down (from elevator trailing edge) 

Vertical tail surfaces 


4° 

20° 


Area (total) 


91 sq ft 


Tin area (two) 


47.8 eq ft 


Span (maximum) 


104.25 in 
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Budder area (two including tabs) 


43.2 sq ft 


Rudder span (maximum) 


101.1 in 


Rudder travel 

Right (from streamline with fin) 


20° 


Left (from streamline with fin) 


20° 


Rudder trim tab area (total) 


3.18 sq ft 


Rudder trim tab span 


36.4 in 


Rudder trim tab mean chord 


6.3 in 


Rudder trim tab travel 

Right (from rudder trailing edge) 


< 


Left (from rudder trailing edge) 


12 



Tub* L ag# 

Maximum width 4 ft 8*5 in 



Maximum height 
Length (tip of nose to tip 
Ingine 

Humber 

Type 

Designation 
lumber of cylinder# 

Gear ratio 



7 ft 4 in 

of tail) 63 ft 5.75 in 
2 

Vright Cyclone 
R-2600-13 or -29 

i 

14 - double row 
16:9 



Supercharger gear ratio 
Low 
High 



7*06il 

10.06:1 



Power settings: 
Condition 
Take-off 
HHP 
Cruise 
Approach 
Tor teste 



rpm Manifold pressure (In 

2600 46 

2400 40 



2100 29 

2000 20 

1800 24 



Propellers 

Manufacturer 

Type 



Hamilton Standard 
Controllable pitch 



— 40 “ 



Number of 'blades 



3 



Diameter 



12 ft 7 in 



Blade type 



6359 A- 18 



Control type 



Hydromat ic full 
feathering 



Huh type 



23150473 



Pitch setting 
Low (fine! 


22° 






High (coaree) 


90° 






Landing gear (estimated data) 








Main gear assembly 


• 


(each) 




Weight of retractable portion 


500 lb 




Radius of gyration of retractable 








portion about retraction pirot 


48 in 






Nose gear assembly 








Weight of retractable portion 


200 lb 






Radius of gyration of retractable 








portion about retraction pivot 


46 in 






Weight 


25,637 


lb 




Moments of inertia (estimated) 










63,000 


slug 


ft 


X ZZ 


120,000 


slug 


ft 


ha 


-1930 


01 

& 


ft 
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FIG A - 1 TEST AIRPLANE B'25J N0 328 




FIG. A - 2. TEST VEHICLE 
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JLPPBTDIX B 



USTHUMETTATION APTB CAL IB EAT I OH 

A. IH STBUMEN TAT I OH 

Ab noted in the Introduction, finding the performance function by means 
of a pulse input haa the advantages of keeping the response within the 
linearity limits of deviation from the equilibrium flight conditions, and of 
requiring a minimum of instrumentation. The devices required for applying 
the input pulses were designed to accomplieh the following functions: 

1. Allow rotation of the control surfaces between any desired limits, 
with positive stops at the ends of the arc of travel. 

2. Allow variations in the size of the pulse. 

3. Allow the pulse to be applied in both directions} that is, allow 
changes in the position of the mechanism relative to the pilot’s 
controls without altering the pulse size. This function also 
allows the pilot to hold two controls rigidly against stops while 
the aircraft is responding to a pulse imparted by the third control. 

4. Allow quick disconnect so as to make it possible to remove all 
impediments from full use of the pilot’s controls in case of 
emergency. 

The mechanisms employed were chosen from many possible designs as the 
simplest that could perform the above functions and still be fabricated with 
a minimum of machining. The aileron pulsing mechanism, shown in Figures B-l, 
B-2, and B-3, consists of a plate mounted on the right control wheel which is 
slotted to receive two pivoted bars which are set to give the desired size of 
pulse. A slotted length of duralumin angle, hinged at one edge, was mounted 
on a shelf fastened to the top of the control column. The slot allowed a 
single stop to be positioned so that the chosen pulse oould be initiated at 
any trim condition by moving the single stop against either one of the stope 
on the wheel plate. The hinged edge allows the stop to be thrown out of 
position in order to clear the aileron control in case of emergency in flight. 
Approximately 45 degrees of wheel rotation were used in applying the desired 
aileron pulse. 
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Ths mechanism used for applying rudder pulsee le shown In Figures 
B— 4 and B— 5. Choice of the design selected was dictated hy the limited 
installation space and by the requirement thAt the pulse setting and quick 
disconnect units must he easily accessible to the co-pilot in flight. The 
mechanism in its final conf igurat ion consists of a long rod, one end of which 
i 8 attached to the right rudder pedal. The other end of the rod slides in a 
block which is piroted at the center, and mounted on a plate fixed to the 
aircraft structure. The block receiving the sliding rod is cut away at the 
top so that a stop pin welded to the sliding rod may move in the slot thus 
provided. The cover of the block is made up of a frame hinged at the inboard 
edge and containing two stop bars connected by a positioning screw to provide 
adjustment of pulse size. In addition, both bars, holding the set pulse size, 
may be moved by means of a second screw to permit engaging the pin on the 
sliding rod with either of the stop bars when in the trim condition. The 
hinge on the inboard side allows the block cover to be moved clear of the rod 
stop, thus quickly and positively freeing the rudder controls. Figures B-4 
and B-5 show the cover block clear of the rod stop in the stowed condition. 
Approximately O.Z in. of linear throw of the rudder pedals provided the 
desired rudder pulse size. 

The microsyn pickoff units installed for measurement of control surface 
rotation were of the variable transformer type designed by Dr. B.I. Mueller 
of the M.I.T. Instrumentation Laboratory. The installation of these units 
in the right and left wing are shown in Figures B-6 and B-7 respectively. 

The mounting positions and linkage lengths are identical in each wing. The 
final determination of linkage sensitivity, that is, piokoff rotation per 
control surface rotation, was governed by accessibility in the wing or rudder 
cavity, linearity limits of the pickoff (approximately 150 mHH-radlans ) ,and the 
maximum control surface motion anticipated in obtaining the desired magnitude 
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of response. Shielded pickoff lends were used throughout, nnd , ns shown 
in Figs. B-6 nnd B-7 , safety brackets were attached to the microsyn mountings 
to prevent any possible damage or Jamming of controls in case the linkage 
worked free at the pickoff rotor. 

The rudder pickoff installation is shown in Fig. B-8 . The linkage is 

i 

similar to those used in the wings, the bell crank at the rudder control 
serving as the initial point of motion measurement. 

The rate gyros employed were elevation and deflection gyros developed 
by the Instrumentation Laboratory, M.I.T. The characteristics of these 
gyros were modified to suit the measurement requirements of the project by 
changing the damping, and were installed as roll and yaw rate gyros 
respectively. 

Laboratory calibration of the gyros prior to installation in the aircraft 
was necessary to determine the excitation currents desired for the elastic 
restraints and the signal generators (pick-offs), the temperature required 
for the damning fluid, and the friction torque level. 

The friction torque level was first determined to Judge whether the gyros 
would provide a satisfactorily sensitive response to allow measurements of 
the desired accuracy and precision. The final results are shown in Figures 9 
and 10 for the roll and yaw rate gyros, respectively. The maximum friction 
torque level for the roll gyro was 50 dy cm, and for the yaw gyro was 140 
dy cme , and the unbalance was reduced to a torque level below that of the 
friction level. The maximum friction torque level for th# roll rate gyro 
represents an angular velocity input of 0.25 milli-radians/sec and for the yaw 
rate gyro, the maximum friction torque represents an angular input velocity of 
0.7 milli-radians/sec. This friction level is considered satisfactory since 
the oscillograph record ordinates can be measured only to a maximum, accuracy of 
.01 inch, the latter figure representing approximately l/2 milli-radians/sec in 
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most of the record® . 



The rate gyro equations of motion were investigated to make certain 
that the desired characteristic could he obtained. The general equation 
representing both gyros is: 

• • • 

W ( ) H8 (llnk)(M M ) ' I (c)(eff) A e * c e A c * 8 (.r)U M ) k c 
g c c c 



where the subscript c refers the indicated quantity to the computer shaft, 
and the subscript g refers the indicated quantity to the gimbal shaft. 

This equation may be reduced to the general form of the second-order 
differential equation yielding 



( ) 



8 



(g)(VA) _ (re) 2 ] + j 2(DE)(FH) 



where 



(DE) 



(FB) 



21 ( V 
c(eff) n 



W 2 = 

\ • n 



(er) (AM) 
X c(eff) 



To determine the gyro undamped natural frequency, the effective moment 
of inertia at the computer shaft was measured using low elastic restraint 
currents and no damping. The period was recorded on a brush type oscillograph, 
and the moment of inertia was then determined from the following equation: 



(eff) 



'(er) (AM) _2 
— T ■ 264 dyne— 



csr-sec 



4ir 



Bechecks at different values of elastic restraint verified the above 
figure. To establish the amount of elastic restraint desired on the gyros, 
maximum rates of roll and yaw for the magnitude of control surface deflection 
inputs anticipated were determined from B-26 roll rate data (~~ ^ ) and from 
Cornell Laboratory Eeports covering fate responses to step inputs of control 
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deflection (Refs. 4 and 5). The mailman rates determined were 0.4 rad/eec 
In roll and 0.5 rad/eec in yaw. The elastic restraint unite were then 
calibrated for both gyros, as shown in Tigs. B-ll and B-12, using various 
values of elastic restraint current. The maximum rate tabulated for each 
sensitivity indicates the rate input for the particular sensitivity which 
will rotate the gyro computer shaft 150 milliradians . At deflections 
exceeding this value, the signal generators become nonlinear. The overall 
sensitivities of the elastic restraints, considering both excitation current 
and angular displacement as inputs are: 



S ( roller) (A, i 2 ; M) 



» 1820 



dyne cm 
(ma) 2 rad 



®(yaw) (or) (A, i\ M) 



5!2 izas-aa- 

(ma) rad 



Selection of the desired undamped natural frequency depends on the 
desired frequency response of the gyro. The rate gyros nust accurately 
reproduce aircraft response below about 10 rads/sec, offering negligible 
or minimum attenuation and phase shift, while at the resonant structural 
vibration frequencies of the airplane, the attenuation must be a maximum. 
Previous instrumentation work on the B-25 has indicated that the most 
pronounced resonant frequencies were found at 16 and 30 cps . To fulfill 
the response requirements outlined above, the following acceptable limits 
for the gyro characteristics were determined: 

Undamped Natural Frequency, n^ » 5 - 8 cps 
Damping Ratio, (DR) * 0.6 - 1.0 

With the above ranges, elastic restraint currents were selected to satisfy 
both undamped natural frequency and maximum rate input requirements. 

A damping fluid was then selected to fit the above requirements of DR. 
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figure B-13 »hove the viscosity vs temperature chamcterietice of damping 
fluid DCC 350, which was the fluid used in both gyro*. DCC 350 offered the 
amalle«t gradient of viscosity to temperature of any of the fluids available 
having the viscosity range desired. 

The following chracteristics , coefficients, and sensitivities describe 
the gyros as used in the final test flights: 





Roll rate gyro 


Taw rate gyro 


x (e)(.ff) 


2 

264 dy cm sec 


2 

264 dy cm sec 




dy cm 


^ dy cm 


c (c) 


17200 nA/.lc 


2400 °Sd/..e 


Damper temperature 


160 ° r 


177° f 


S (.r)(AM) 


3.8 I 10 s 

rad 


4.66 X 10 6 

rad 




38 rad/ sec 


42 rad/sec 


n 

n 


6.05 cps 


6.7 cps 


(DR) 


.857 


1.08 



The gyros were mounted in a dustproof container and secured in the air- 
plane on an adjustable mount which allowed aligning the gyros with the 

horizontal. Since the B-25 trims at 175 mph in a nose high attitude with 

o 

its horizontal reference line 3.6 above the horizontal, the gyro container 
must be tilted accordingly. The adjustable mount is shown in Fig. B-14, 
and the gyro container in position on the mount is shown in fig. B-15. The 
gyro container utilized three thermostatically controlled heaters to maintain 
a fixed ambient temperature for the gyros. 

figure B-16 shows the amplifier panel, four channels of amplification 
were used to amplify the piclc-off voltages representing total aileron deflection 
rudder deflection, angular velocity in yaw, and angular velocity in roll, 
lach amplifier channel performed three functions (l) amplified the AC picfc-off 
signal, (2) demodulated the AC signal, and (3) amplified the resulting DC 



signal. The output of each amplifier channel is controlled by 12 attenuation 
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setting* in order to adjust oscillograph current to give the desired 
oscillograph displacement. 

The oscillograph used is shown in Fig. B-17. It is a 12-channel 
Consolidated Recording Oscillograph. The dynamics of the oscillograph 
galvanometer units have a negligible distortion effect in reproducing the 
inputs introduced hy the amplifier outputs. The natural frequency of the 
units used is 120 cps and the DR is 0.7 approx. 

The stiffness current regulator shown mounted in position in Fig.Iarlfl 
is a closed-loop feedback regulator functioning to hold the stiffness current 
at a set current within 1 percent variation. Figures B-15 and B-17 show the 
exposed regulator chassis, the chassis having been removed from its case. 

To integrate the added instrumentation with the instrumentation 
previously existing in the aircraft, an additional Junction box was 
installed. It’s location i6 noted in Fig. B-17. 

The three-phase, 400-cycle thyratron controlled inverter initially 
installed was found to give unsatisfactory regulation for the needs of 
this project. Consequently, it was replaced by a USAT PO-16 single-phase 
inverter, which was found to give satisfactory regulation. This change 
necessitated use of a phase splitter and stepdown transformer to obtain 
the 3 phase, 28v gyro wheel supply. A block diagram of the overall 
instrumentation system is shown in Figure B-26 . 

As shown in Fig. B-26 the output of each pickoff was loaded with a 
1000-ohm resistance to reduce the noise level and the null voltage. A 
variable resistance and a variable inductance were placed in series with 
the pickoff excitation coils (P.0, coils in series). The resistance was 
used to set the desired excitation current, and the inductance was used to 
bring the pickoff outputs into proper phase. 
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B . CALIBRATION 



The objective of the measuring system calibration for this project 
was to establish overall sensitivities, expressing Inputs to the measuring 
instruments in terms of output linear displacements at the oscillographs. 

The sensitivities are expressed in such a manner that multiplication of any 
displacement on the oscillograph record by the appropriate sensitivity yields 
the measured quantity in the desired units. These sensitivities are: 



(all) (x a ) 




deg/ in 


osc a 






(rud) (I a ) 




deg/ in 


osc r 




deg/ sec 


(roll gyro)(X W_) 

QIC A 


in 




v ) 

V 


deg/ sec 


(yaw gyro) (X 

0 8C 


in 



Tor purposes of simplicity and system flexibility, the overall sensitivities 
listed above were broken into two component sensitivities, one representing 
the instrument sensitivity and one representing the indicating system 
sensitivity. The instrument sensitivity is a fixed quantity determined in 
the case of the control surface rotations by the linkage sensitivities aDd 
the signal generator excitation current, and in the case of the rate gyros 
by the elastic restraint current and the signal generator excitation current , 

The indicating system sensitivity may be varied at will by changing the 
attenuation setting at the amplifier. The amplifier sensitivity allows the 
desired size of recorded input pulses and output angular velocities to be 
selected. Selection of the proper size of recorded displacements is governed 
on one hand by the linearity limits of the oscillograph trace and on the 
other hand by the size record required for ease of analysis. 

In order to obviate the necessity for calibration of each amplifier at 
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•ach attenuation setting, a ’’master'' pickoff and a calibration voltage 

source, the latter providing 12 voltage etope, were employed. The master 

pickoff was calibrated (Tigure B-18) to determine the limite of linearity. 

The sensitivity determined was 3/ w. v e 4 . with 300 

(■poHA e ) division 
mpo mpo 

master -p.o. divisions as the limit of the linearity range. The master 
pickoff was then used to ’’buck-out’’ or neutralise the voltage produced 
by the instrument plokoffs. Thus, instrument sensitivities were obtained 
both in terms of voltage output for input rotation or angular velocity and 
in terms of master pickoff output angle for input rotation or angular velocity. 

The master piokoff is shown in Tig. B-19. The master pickoff could be 
plugged into any of the pickoff outputs by means of bayonet Jack receptiole* 
on the amplifier panel. 1* shown in Fig. B-26, the particular pickoff out- 
put could then be switched from the amplifiers to the master pickoff. 

The procedure for obtaining the instrument sensitivities was straight- 
forward. Tor the ailerons, a propeller bubble protractor, measuring to O.i 
degree accuracy was used to measure the angle of both ailerons with respect 
to the horizontal, and the total aileron differential angle thus obtained. 

Tor each measured ft , the corresponding voltage was recorded. This voltage 
was then nulled by the MPO and the corresponding angle of the MPO recorded. 

The voltage representing total aileron angle was achieved by adding the 
voltages of the two aileron pickoff outputs in series. 

Tor calibration of the rudder deflection to MPO angle, the same system 
was used as outlined for the ailerons with the exception of the means of 
measuring rudder deflection. Figure B— 20 shows the method employed. The 
center of rotation of the rudder was determined by means of a plumb bob. 

A pleoe of mirror was then attached to the rudder at the center of rotation, 
and the plumb bob lined up with the center of a sector of a circle, the 
latter being graduated in mils. Measurements were then read by means of a 



- 53 - 



telescope equipped with a hair line. 



Figure B-21 shows the set-up used to calibrate the gyros. On a small 
turntable, having a selection of 4 speeds, an inclined plate was mounted to 
receive the gyro container. The can was so aligned that the spin axis of 
the roll rate gyro was rotated a measured angle about the gimbal axis. With 
the two gyros properly aligned within the can, this also rotated the input 
axis of the yaw gyro about its spin axis by the same measured angle. As a 
consequence, the input to each gyro was calculated in the following manner: 
Roll gyro sin 6 

Taw gyro * z = * TT cos € 



where is the angular velocity 
of the turntable, and 9 is the measured angle of tilt. 

The calibration curves for all four measuring instruments showing the 
measured quantity vs MPO divisions (each div. = approx. 2 mils) and, for 
the gyros, the measured quantity vs pickoff output voltage are shown in 
Figs. 2-22 through £-25, 

The instrument sensitivities thus determined are 
[aUeron] S (all)( , # ^ S (Mp0)( ^ %) - S ( all ) (X) (A^ «.) * °' 0281 Ilf 



[rudder] 



[roll gyro] 



[yaw gyro] 



S (rud)( V J r ) S (KPO)(i mo .„> = SfrudHtOU^ «,) 
S (rg)(« o W x ) S (MPO)(A Mpo « o ) = S (rg)(M)(A >!K1 Wj) ‘ 

S (y«)(«o V Z ) S(MPO)( Vo •„> = S (« )(#,( Vo V ‘ 



0.00846 



deg 

div 



0.994 



mils/ sec 
div MPO 



1.228 



mils/ sec 
div MPO 



The indicating system sensitivities were determined on each flight 
whenever data was recorded. The indicating system sensitivity is made up 
of the amplifier sensitivity, oscillograph sensitivity and master pickoff 
sensitivity, that is: 
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S (ind aye) (l Q9C . A^pg) " S (MPO)(e A^pg) S (ampXi,e) S (osc)(l ogc 0 

Upon the completion of each series of flight rune the measuring 
instruments were switched away from the inputs to the amplifiers, and a 
12-step voltage source switched onto the amplifier inputs (shown in Fig.B-19). 
As the voltage was varied through the 12 steps, a record was made on the 
oscillograph. These same 12 voltage steps were then applied to the master 
pickoff, the latter instrument being rotated until each calibration voltage 
was ndlled and the corresponding output angle of the master pickoff was 
recorded. Thus increments of oscillograph displacement corresponding to in- 
crements of MPO angle were available for each test flight, and an overall 
indicating system sensitivity was available without knowledge of the 
oscillograph and amplifier component sensitivities. That is, 

A A MPO 

(lnd .y.)(X otc A^) ' A X 0 , c 

The overall system sensitivity for each measuring channel was thus 
available by combining instrument and indicating system sensitivities. 
Therefore, the general expression of the sensitivity ie 

8 (l ehanH* oee q) * S ( XMXA^, q) S ( Xlnd .y=)d oec A^) 

To determine the system errors and uncertainties, static teste were 
made on the ground. A known value of eaoh quantity to be measured was 
applied to the respective measuring instruments, and an oscillograph 
record made of the outputs. The calibration source voltages were then 
recorded at the same attenuation settings, and a record of master pickoff 
angles reoorded for each calibration voltage. In this manner the overall 
sensitivity could be compared with the product of the instrument sensitivity 
and the indicating system sensitivity. 
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The comparison was found to be; 
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The linearity limits and of the measuring system components are as 



follows: 



Component 
Maeber P.0. 

Aileron Defl. 
Rudder Defl. 

Roll gyro 
Taw gyro 

Oscillograph defl. 



Linearity Range 
<500 divs ( 150 mils) 
30° total deflection 
^7° rudder deflection 
40.3 rad/sec 
40.35 rad/ sec 
42.7 in (approx.) 



The largest source of error in the measuring system lies in the 
insensitivity of the voltmeter used as an indicator when nulling a voltage 
with the master p.o. Ifforts to reproduce results have produced variations 
as high as 4 divisions of the MPO, with an average error of 1 to 1 l/2 
divisions. This represents an error, or null voltage, of 4 to 6 mv in matching 
voltages. 

A change in the instrumentation was effected Just prior to the final 
test flight (See Flight Record - Appendix D) . The gyro container was 
secured rigidly to the aircraft structure during flight in order to eliminate 
recording vibrations set ip by the gyro container mount and the Lord shock 
mounts. This modification resulted in removing an undesirable 10 cps 
transient from the records which occurred whenever the aircraft was respond- 
ing to input pulses. This change made it possible to make the shock mounts 
ineffective during flight, and yet to free them when desired in order to 
protect the gyros from landing shock. 
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MECHANISM IN STOWED CONDITION 
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rn^TROi DF'/irp Ann v v ai. ERON PJLSF - RIGHT SIDE VIEW 
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CONTROL DEVICE FOR APPLYING RUDDER PULSE 









CONTROL DEVICE FOR APPLYING RUDDER PULSE - TOP VIEW 
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FIG.B-9 UNCERTAINTY MOMENT VS COMPUTER SHAFT ROTATION FOR 

ROLL RATE GYRO. 
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FIG. B-IO. UNCERTAINTY MOMENT VERSUS COMPUTER SHAFT DEFLECTION FOR YAW GYRO 
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APPLIED MASS IN GRAMS 
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ADJJSTABLE MOUNT FOR GYRO CONTAINER - 
VIEW LOOKING FORWARD 

re 
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15 LOCATION OF GYROS - LOOKING FORWARD IN AFTER COMPARTMENT 




AMPLIFIER panel 
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MASTER PICK OFF 



MASTER PICKOFF 
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-19 MASTER PICKOFF AND CALIBRATING VOLTAGE SOURCE - LOCATED 

BELOW AMPLIFIER PANEL 





FIG B- 21 SETUP USED FOR G^RO STATIC CALIBRATION 
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APPENDIX C 



ANALYSIS OT A DIVERGENT TRANSIENT RESPONSE 
Reference* (2) and (3) point out that the approximation of the Fourier 
Transform is applicable only to stable system*. This can be easily *een from 
a consideration of the transient response of a system described by a 
characteristic equation haring a positive root. For such a response, the sum 



I 

n 



1 , 2 ,. 






will not converge. This condition is found in the lateral motion of as my 
airplanes, where the positive root is associated with the motion called 
spiral divergence ♦ 

Since the B-26J used in the present investigation did not exhibit this 

type of response, a simple dynamic system illustrating divergent motion is 

» 

analyzed in this appendix, using the finite pulse method. 

Consider the system schematically pictured in Fig. C-l . The mass, m, 
in this figure is assumed to be supported in the gravity field by a rigid, 
massless, rod. The Laplace Transform of the characteristic equation for the 
small angles of motion considered in Fig. C-l, (Ref. 8), is 
[m-^ s^ ♦ C ■£ s - ng] ® ° 0 (C-l) 

Assume that a forcing function is applied consisting of a unit triangle 
of force, UT(t - At) which is described by 
T ® 0 0 > t 



F =* — t 0 ^ t 6. dr 

At 

F * 2 - t*- t At = t = 2 At 

At 

F - 0 t > 2 At 

At * 0.1 second 
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The transformed equation of motion becomes 

[ml s 2 + C _1 ■ - mg]®(s)b (LT)[UT(t - At)] 

where 

(LT)[UT(t - At)]«= S±- ~ J 

At e 



Solving for ® (■) , 
® (a) = 



/ _ — At fc \ 2 

Al ~ e ) 



(At a^)(ml)(e^ + — 




(C-2) 



(03) 



In this investigation of the B-25J lateral motion, the response which 
was measured was the angular rate in roll or yaw. For this reason, the 
time rate of change of the angular position® will be considered as the output 
in this example. The rate of change of this output,® , is given by 



®<e> = 



(1 - e" 6 ™) 2 



At mi i(i + ijj— s - J-) 
Factoring the characteristic equation, 
C 



(04) 



2d -' 



,2 

4m* 



C — + ■&— 
2 j 



Let 



2m 



♦ <i 



c ^ £^_ 
I? i 




Bote that a is always positive and that a and b are positive, real 
numbers . 

By appropriate substitution of a and b, the equation of motion, eq . (C-4) 
can be written, 

®(s) = ^ (C-S) 

At mls(8-a)(s+b) 



- 84 - 



The performance function (P F) 



is defined as 



[F.0 ] 



(PF) 



[T,0] 




• - J» 



Then* 



(PF) 



[ t . i ] 



1 Je> 

ml ( - a) (Jut) ♦ b) 



(C-6) 



A value of 0.1 was assigned to a, and a value of 2.0 to b. The value 
of 0.1 for a is such that the amplitude of the divergence will double in 6.93 
second*. Some airplanes have this much divergence at high angle* of attack. 
The valu* of 2.0 for b i* somewhat small compared with that for the highly 
damped notion of nost airplane*. The major consideration in the choice of 
these roots was to choose values that would reduce the calculation* as nuch 
as practicable, and yet that would retain some eemblanoe to those encountered 
in airplane lateral performance functions. 

These choices were consistent with the choice of m and i as 1 and 2 
respectively . The units are those of any consistent set of physical units. 
Substituting these numbers and letting 
s ■ J® 

the performance function becomes 

1 J© 

(Pr) [T. 0 ] “ T (job - O.lHjui +" 2) (07) 

This performance function is plotted in Fig. C-5. 

Returning to eq . (C-5) , by the inverse Laplace transformation the time 
response was obtained and plotted in Fig. C-2. 

In Fig. C-3, the portion of the curve of Fig. C-2 beyond 6 seconds was 
plotted on semi-log paper, thus enabling determination of the divergent root. 
The straight line was extended to eero time and the intercept found. Having 
thus found an expression for the divergent component of the response, it was 
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subtracted from the total response leaving the remainder shown in Fig. C-4. 
The integral of this remainder is convergent and its approximate Fourier 
Transform is 

(rr)0. Ax) = (rT)|oT(t3 £ ®(t 

(r,B) n = 1 , 2 ... n 

Since the system is linear, the performance function determined from 
the response curve can be expressed as the sum of components, i.e., 



(W) 



[T, ® ] 



(LT) ^o (t) ]|dlverxe nt) 



(FT)Ca (t)3 



~mr< 



(rsm) 



in 



(TT)fq in (0' 



(rt)t(RT)[q (t )], , 

o o J (rem) fn _ n \ 

(M|V(t)] ( °- 0) 



in 

The last two terms on the right side of eq. (C-8) were determined by- 
approximating the Fourier Transform by triangles. The first term can be 
determined since the analytical expression for q^t)^^ is known, and 



since 



CLT)[q ln (t)] = (TT)[q ln (t)] 



This analytical expression for the divergent component of the response as 
determined from the time response is 



UT>[q o <Od lT ] 



0.0231 
s - 0.1 ’ 



and 

. 0.0231 

[F, ® ](diT) (j« - o.l)(FT)[OT(t)]e"' ,4T ” 

From this, the amplitude ratio and phase angle can be calculated for any 
given w . 

The desired performance ‘function can be obtained by performing the 
addition of complex quantities indicated in eq. (C — 8 ) • 

Figure C-5 shows the comparison of the original performance function 



with the performance function obtained by analysis of the time response. The 
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acouraoy with which the original performance function wae duplicated shows 
that this ie a useable method of determining a performance functibn from a 



divergent time response 



VISCOUS 

FRICTION 

COEFFICIENT 




FIG. C-l SCHEMATIC DIAGRAM OF SIMPLE DYNAMICAL SYSTEM 
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SECOND 




RATE OF CHANGE OF ANGLE , © , IN RAD/SEC. 
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FIG. C-5. COMPARISON OF THEORETICAL AND MEASURED PF 



APPENDIX D 

DATA REDUCTION AND EVALUATION 

The method used to obtain a deeired performance function, knowing the 
time response of a system to a known input pulse, is outlined in Chapter I. 

In that chapter the method of approximating the mathematical expressions 
graphically (within the deeired degree of accuracy) is cowered. Reference 3 
shows the method applied to a critically damped second order system, and out- 
lines a method of attack to be used when more lightly damped oscillations 
are present in the time response. The process as abore outlined has been 
applied to the oscillograph test records obtained in this lure stigat ion. 

The oscillograph records provided the following information: 



Input Pulses 

Aileron deflection, (+ defl) 
Aileron deflection, 6 (“ defl) 

A 

Rudder deflection * (+ defl) 

Rudder deflection, fl (“ defl) 

r 



Output Response 
Angular rel. in yaw and roll 
Angular rel. in yaw and roll 
Angular rel. in yaw and roll 
Angular rel. in yaw and roll 



Sign conventions used in designating pulses correspond to NACA usage , i.e., 

(«*■) Aileron defl. - right aileron depressed 

(♦) Rudder defl. - rudder trailing edge moved to the left. 

° Reduction of the above listed recorded data would yield the following 

aircraft performance functions: 



(1) 


C^(A)(. . 

A 


. * x ) 


(2) 


CPT] U)<V w z > 


(3) 


[P, 1(A)(V 


V 


(4) 


tP, lu>(« r . w z ) 
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The time available for redaction of data was not sufficient to allow 



all these performance functions to be determined. Therefore, since this 
project is primarily concerned with the application of the pulse technique 
as a means of determining the performance function of any dynamical system, 



and since knowledge of the actual dynamics of the test vehicle is of 
secondary importance, the records were analysed to achiere the following 
objectives: 

(1) To determine the degree of accuracy or uncertainty in repeating 
results. 

(2) To compare the results of determining the performance function by 
slightly variant graphical means. 

(3) To compare the performance function as obtained by the pulse 
technique with the calculated performance function as determined from the 
aerodynamic derivatives obtained from wind tunnel teste. 

(4) To reveal any inherent uncertainties or limitations in using the 
pulse technique to determine performance functions. 

To obtain the data necessary to investigate the above objectives, the 
following analyses were actually made: 



TABU 1 



Bun Ho* 



Method 



Direction of 
Pulse Input 



4665 



4661 



4662 



4662 



4666 




CW] t‘a V 



Graphical 




Graphical 



Graphical 



Graphical 



(♦) defl. 
(-) defl. 
(-) defl. 
(+) defl. 
(+) defl. 



4665 



Numerical 



(■*■) defl. 
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In determining the performance function graphically, three different 
methods were utilized to determine the amplitude ratio and phase angle at 
each selected forcing frequency. The first method utilized the general 
method as outlined in Ref. 3. that is, the graphical summation of input and 
output vectors as indicated in the following expression: 



[FF] = 



£q/ ,\(t )e 

(out) d 



-JnwJW 

0 [FT][UT(t)] 



-Jnw-AT . 

Iq (in) (T n )e [FT][UT(t)] 1 



(D-l) 



n = l|2|3f4|5| 

Figure D-l shows a sample plot of this method for w^ » 1.25 rad/sec using 

run number 4665, for determining the performance function [PF]^^^ W )* 

a X 

this plot, the effect of the lightly damped sinusoid on the phase angle and 
amplitude ratio is apparent. Each cusp represents one half cycle of 
oscillation. In order for the amplitude to be attenuated to 0.05 of the 
original amplitude, the plot must be continued until a total of approximately 
6.5 cusps have been plotted. The majority of the oscillograph records are 
not long enough to allow more than three or four cusps to be plotted, not is it 
practicable to obtain transient responses of this length with a usable 
uncertainty level. Thus the AR or PA, or both AR and PA are subject to 
error due to deleting the vector summation representing several additional 
half cycles of the oscillatory mode. Furthermore, the desired matching of 
the curve required that the response be divided into 0.1 sec increments, 
requiring that as high as 84 vectors be plotted (8.4 sec) at each forcing 
frequency — a tedious Job resulting in a rather high degree of uncertainty. 

In order to simplify the graphical method by requiring fewer vectors 
to be plotted to determine the AR and PA at each frequency, and to obviate 
the error incurred by omitting two or three half cycles of the oscillatory 
mode, a second graphical plotting system was used. In this system vectors 
were plotted in identically the same manner as previously discussed until 
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the time response was reduced to a remainder which was solely (for practical 
purposee) a pure damped sinusoid. At this point a single vector was added 
to the response vector summation to represent the contribution of the long 
’’tail” of the sinusoid. 

Figure D-2 is a sample plot of this method at the earn forcing frequency 
and run (No. 4665) as used in Fig. D-l. The accuracy of this method depends 
on the accuracy with which the characteristics of the damped sinusoid may 
be determined. The characteristics were determined by the method demonstrated 
in 8 of Ref. 3. Knowing that the oscillatory motions in both yawing and 

rolling angular velocity have the same frequency and damping ratio 
(see eqs. (D-10) and (D-ll) several records of both motions were used to 
measure the period and peak amplitude ratios, and average values for the 
parameters, DR and W^, were thus obtained. The logarithmic decrement curves 
shown in Fig. 7 of Ref. 3 were used in determining DR. The records agreed 
quite closely in the period of the oscillation, but determination of DR 
required the average of several records, giving the following characteristics: 



DR = 0.15 



W = 1.45 rads/sec 
n 

In the case of the rolling angular velocity response, it was found 
that the response was nearly a pure sinusoid after the third peak overshoot. 
(See Fig. D-5) . As a consequence, the general performance function equation 
was rewritten, designating the time at the second peak overshoot as t^, and 



the ordinate at t as 



[FT] 



P a ’V t (£ 1 (out) (T n )e 



■Jnw Ay 

1 0 [FT][UT(t)] o ; 



( q out^ 



n ~ ~ 1 ,2,.., p 



-jnw^Ai 



1 ) [W]tOT(t)] ] 



[(£ ^(in)^ T n )e 

-j dw/t ) 

[(Iq (t )s r ° [FT][UT(t)l 

on o 



1 , 2 , ... 



i 

b = p,p + 1 ♦ . . 



-.Jnw JW 

[(Iq, (Oe ^ 4 ) [FT][UT(t)].] 

1 n 1 In = 1,2,... (D-2) 
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(zero time lag, t , is assumed, in the preceding equation) 

The second term in the numerator may now he replaced by a single rector 
representing the contribution of the damped sinusoid from t^ to infinity, 



yielding, 



“Jnw At 

2 »( 0 »t> (T n ># ° 



In ^out 



“Jnw^AT 



n ^ 1 1 2 [ « » 



2 






[fi][nr(t )] 1 

-JpVfA T 0 



n = 1,2,.., 



(D-3) 



Z q 



(in) 



(r n ) e 1 [FT][UT(t)] i 



n — 1,2,3,.., 



where and are the time intervals used in the output and input 
respectively . 

Any time later than the time at which the response becomes a pure sinusoid 
may be selected as the point at which the remaining oscillatory motion may be 
represented by the Fourier Transform of a damped sine or cosine wave. However, 
it i« advantageous to use a * 'peak point’’ or a point where the sinusoid 
crosses the axis in order that the phase angle in the sinusoid expression may 
be eliminated. In these two cases the term is reduced to the following Fourier 



Transform expressions: 
(t at peak) 



t fT K(osc) 






“DBw t 

a *. 

e cos w t 

o 



Jp V T . 



1 + / FR \ 2 
V DB ; 



tan 



-1 



,FHx - tan" 1 2(DB)(FB) 
W 1 -IB 2 



- (TB ff * [2(DH)(ra)] : 






(D-4) 



(t at axis) 
P 



[FT]q 



o(osc) 



[W]qp 



-(DB) (w )t 



n 



sin w t 
0 



•Jpw^i 



q w 
o o 



C 



-1 



“tan 2(DB)(FB) 

1 - ft&Z 



^tl “ (FR) 2 ] 2 v (2(DB)(FB) 






> e 



'Jpv f AT 1 



(D“5) 
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The time, t^ should be chosen to fall on one of the ordinates determined, by 

the choice of At^ to avoid disturbing the area aporoximat ion originally 

established by the isosceles triangles. If a peak maximum or minimum is 

selected as t , the area approximation by triangles is interrupted even 

through t may lie on one of the ordinates of the fitted triangle. In this 

case, a right triangle term must be added to the Fourier Transform of the 

output if the resulting output vector i 6 to be exactly correct. However, 

the initial amplitude of the cosine wave starting at t was small enough 

P 

in the runs analyzed in this project to render the right triangle term 
ineffectual in changing the AR or PA of the performance function. The 
method is demonstrated in Fig. D-4. 

Figures D-l , D -2 and D-3 compare the two graphical methods used in 
determining the performance function, £****] ( A ) (5 \f )• D-l shows 

the AR and PA determined by plotting the maximum number of vectors 
available from the oscillograph record. Figure D-2 shows the same AR 
and PA at the same forcing frequency as determined by representing the 
output divided in two parts in the method just described. Figure D -3 
shows a comparison of the AR and PA results determined by each method. 

It was found that the maximum error in the first method, Fig. D-l, appeared 
at the resonant frequency. 

The third method used for determination of the performance function 

from pulse response data is the method described in Ref. 3, Section IV. 

The oscillatory mode was determined as previously described, and 

subtracted from the output response from t = 0 to . As shown in 

Fig. 9 of Ref. 3, the remainder of the output response is no longer 

zero at t = t , and a term must be added to account for the right triangular 
o 

area at the time origin, which is not included in the area approximated by 
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the iiosceles triangles approximating the remainder of the response. The 
resulting expression is that given by eq. 40 of Ref. 3. 




in out 




) ( t ) ^[Vrem) (T o ) ^ rT ][ (UBT)(t) 3 



2 

n = 1,2 



• • • • 




t • • • 



■jnv At 

2 £ q i (T n^* ^ 1 C^T]C (UT) (t)D 

n = 1,2,... (D-6) 



» • • • 



Figures 3-1 and 3-2 show the amplitude ratio and phase angle 
respectively of the performance function for aileron deflection input and 
roll response output. The phase angles shown have "been corrected for gyro 
phase lag as obtained from the Fig. 0-13. Undamped natural frequencies and 
damping ratio for both gyros are found in Appendix B. In figures 3-1 and 
3-2 three experimental performance functions are shown which are the result 
of analyzing graphically three different oscillograph records (runs 4665, 
4661 and 4662) . The uncertainty in reproducing results may be observed 
from these plots. The curve shown represents the average value of the 
three experimental points at each frequency. The greatest deviation from 
the average occurs at w^ * 1.25 where the maximum variation in AR is 5.9 
percent. The average variation is 1.6 percent. As noted in the sample 
oscillograph records, Figs. D-5 and D-6, no two input pulses were identical. 
Consequently, each response is different which necessitated independent 



analysis of each record. The reliability of the pulse input technique in 
reproducing results appears to be excellent, judging from the results of 
this limited investigation. 

From the experience gained in use of the analyzing systems previously 
outlined, the following advantages and disadvantages were noted in applying 
these systems to time records having a long ’’tail*’ due to a lightly damped 
sinusoid . 
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METHODS OF ANALYZING THE RESPONSE 
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The only remaining curre shown in Figs. S-l and 3-2 which has not been 
discussed is the calculated performance function using wind tunnel stability 
derivatives for the B-25. Due to differences between the model used in the 
wind tunnel teets and the actual airplane, these derivatives have been 
modified to more closely represent the true derivatives of the aircraft. 

The derivatives that were not available from wind tunnel tests were estimated 
on the basis of previous engineering experience with comparable aircraft 
configurations. For use as a means of comparison in this project, the 
derivatives were corrected to include the values of moment of inertia about 
the I and Z axes as found for the B-25 used in the testing configuration 
(See Appendix A), 

The derivatives as used to determine a calculated performance function 
are as follows: 
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* -3.175 l/sec 2 
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* +.00648 l/ft-sec 
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The sign convention for the derlratiraa and aquations of motion and 
the symbols representing displacements, moments, velocities, etc., are In 
accordance with NACA designation system. Positive directions along the 
axes, viewed from the origin, are forward, right wing, and down for the 
X, T, and Z axes respectively. 

The following assumptions were made in writing the equation of motion 

1. Damping moments due to control surface motion are negligible. 

2. The wind axes and the aircraft reference axes are considered to 
be superimposed. 

3. No gyroscopic effects were considered. 

The equations of motion are: 

a. S umma tion of forces along Y axis 

nrf ♦ mUr = T + ^ (D“7) 

b. Summation of rolling moments 

%k T ♦ „ «. ik r ♦ t ; + , 






dr ' dp r dr U * ' e«_ a 



c. Summation of yawing moments 

i„ ; . T . *e + m + 

dr 
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*2 ‘ ‘ eV ’ T 55 p ’ 5? r * hz * * H, *« 



(D-8) 



(D-9) 



Dividing each equation by the appropriate moment of Inertia or the 
naee, talcing the Laplace tranef ora and solving for p and r, the following 



performance functions were obtained: 



[»](« v_) 

ft 



■3.175 S (S 2 v ,5863 + 1.6) 

(S 2 v .377 S v 1 .78) (S + 2.9)(3 - .0013) 



where p « W^. 





.184 (S ♦ 3.5l)(8 + 1.097)(S - .837) 

(S 2 + .377 S + 1.78) (S v 2.9)(S - .0013) 



(D-10) 



(*- 11 ) 



where r = W^. 

The calculated performance function [FF]/^ v ) shown in Tigs. 3-1 

a I ; 

and 3-2 have the same general shape as the experimental performance 
function. The misalignment of the resonant frequencies of the two plots 
is indicated in the difference in characteristics of the oscillatory modes 
of the calculated and experimental performance functions. The characteristics 
for the calculated and experimental oscillatory modes are: 



DB ¥ 







n 


Experimental 


.15 


1.45 


Calculated 


.141 


1.36 



Since the oscillatory mode is the predominant mode in establishing 
the peak values of the performance function, the difference in the 
magnitude of the two peaks should represent approximately the same 
percentage difference as indicated by the difference in their damping 
ratios. From Tig. 3-1 the peak of the experimental performance function is 
93.8 percent of the peak value of the calculated performance function, while 
the corresponding ratios of the two DR ie 93 percent. 

The calculated performance function shows a very small divergent root. 
However, it is difficult to state whether the B-25J actually is spirally 
divergent or not. The oscillograph time records do not Indicate a divergence. 
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However, in order to determine whether suoh divergence exists, it would "be 
necessary to take a record of at least 15 seconds to determine the presence 
of this unstable mode. During this period small perturbations in roll or yaw 
due to disturbances of the air mass would probably nullify the accuracy of 
any such record. In any case, for practical considerations, a root as small 
ae the calculated performance function indicates, whether positive or 
negative, would have negligible effect on the performance function. 

The sensitivities were compared using the following assumptions: 

1. The unstable root encountered in the characteristic equation of 
the calculated performance function was considered negligible 
(0.0013 * 0), 



2. The experimental performance function contains no divergent roots. 

(An assumption implicit in positioning the base line.) 

Using the first assumption, the sensitivity for the calculated response was 

determined by rewriting the performance function as a sensitivity multiplied 

by a frequency function. The same assumption applied to the experimental 

curve gives the ratio of the area under the response curve in the time 

domain to the area under the input curve in the time domain as the sensitivity 

of the experimental performance function. 

To determine the area contributed by the oscillatory mode of the 

reeponse extending beyond the limits of the oscillograph record, the area 

was measured out to t * t , (t again representing the point at which the 

P P 

response become* a pure damped sinusoid) and the curve beyond time t^ was 
integrated analytically and added to the area from t =0 to t^. The 



expression thus becomes 



[I q (t ) At ] 
on o 



S (A)(6 W ) 

a X 



a X experimental 



t = °° -PB V t 

♦ / %(t) • n 

t = p o 

n * -Ll.3 1 »-• »p - 

vV &T i ] 

n = 1,2, ... 



COS V t 
0 
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The sensitivities thus determined are: 



^'•aVteale) 



= -.992 1/ sec 



,<A)(8 a H lW r ) 



* -.945 l/sec (Bun Ho. 4666 Input • ♦ 6 ) 

A 

■ -1.080 l/sec (Bun No. 4666 Input = + 6^) 

= -.985 l/sec (Bun No. 4661 Input = - fi ) 

a 

= -.990 l/sec (Bun No. 4662 Input = “ 6 ) 

a 

An additional sensitivity for comparison with the above may be obtained from 

Pb 

flight test data for the B“25, utilizing curves of vs 6^ (Bef. 4). From 
these curves, the sensitivity is, 

S(A)(J a V W (test data) = '- 9 ° 

The comparison of experimental and calculated sensitivities are in 
close agreement, as is the check with flight test roll rate data (Bef. 4) . 

The performance function for yawing velocity response to aileron input 
is shown in Figs. 3-3 and 3-4. The yawing velocity response as shown in 
Fig. D-5 appears to be an excellent example to illustrate the limitations 
of the pulse input technique. As noted in Fig. D-5, the response is 
essentially oscillatory, the first lobe (positive) indicating the adverse 
yaw effect in which side slip has taken place to give the high righting 
moment responsible for the second lobe (negative) being of greater amplitude 
than the first lobe. The pattern of this response is noted in detail since 
its effect in the graphical analysis is important. 

The type of output response pattern found in yawing velocity for aileron 
input shows the same summation pattern at each of the low frequencies. Com- 
paring Fig. D-7 with Fig. I>-2, it is seen that the end point of each 
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summation is well away from the origin, and that small errors made in 
plotting will offer very little error to either the Ail or the PA of the 
total output vector. Also, the comparison at these low frequencies with 
the AH and PA of the calculated performance function is satisfactory both 
in yaw and in roll. (See Fig. 3-1 and 3-2 for roll, and Fig. 3-3 and 3-4 
for yaw.) However, as the frequency increases, and W^At increases corres- 
pondingly, the vector summation of the yaw vectors begin to ’ ’wind up,’’ 
the end points never getting very far from the origin. Note the difference 
in the summation as shown for = 1.5 rad/sec in Fig. D-7 and the summation 
as shown for = 5.0 rad/sec in Fig. D-8 . This clearly indicates that any 
uncertainty entering the summation, regardless of its nature, will impart a 
much higher uncertainty effect to the PA of the output vector (end-point to 
origin) than in the case where the end-point is well removed from the origin. 
The uncertainty of the results thus are increased by two factors as the 
frequency increases, 

1. The large relative size of the last vector in the summation in 
comparison with the resultant output vector, and 

2. The accumulated errors from measuring ordinates from the record. 

Consequently, each system of analysis as covered in Table 2 must be 

analyzed and tested to determine the relative uncertainty in approximating 
the Fourier Transform under conditions where the response is composed in the 
main of an oscillatory motion. The general effort is to establish the system 
of analysis that will produce the fewest inaccuracies in approximation and the 
fewest uncertainties in plotting. 

Investigation of the method listed in Table 2 as (2), the one Used 
most extensively in determining the response vectors for the performance 
function in roll for an aileron input, certain inherent uncertainties are 
noted. Again, it is emphasized that these uncertainties become increasingly 
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important as the response becomes more predominately oscillatory. These 
uncertainties affecting the magnitude and phase angle of the output rector 
are listed in order of the importance of the error they present. 

1. Uncertainty in determining the time at which a mode has reached its 

peak value on the time response, or the time at which the output time response 

crosses the base line. That is, correct determination of the location ef t . 

P 

Figures D-5 and D-6, actual oscillograph records, show the cause for this 
uncertainty. The records are of too short duration to show more than a maxi- 
mum of one full cycle of the pure damped sine wave. Therefore, although the 
DR and W n of the oscillating mode may be well established by comparison of 
many records, the low flat-topped lobes and the low slope of the curve, plus 
the fact that the noise to signal level is becoming quite high at low 
amplitudes, all tend to make the selection of the exact point at which the 
curve peaks or crosses the base line somewhat uncertain within about +0.1 
second. The effect of such an error can be observed best at high frequency. 

If W = 6 rads/sec and&T = 0.1, and the final output vector is approximately 
the length of the vector representing the oscillatory mode, an error as high 
as 30 degrees in PA may exist. As the end point moves further from the 
origin, this possible maximum error reduces accordingly. 

2. Uncertainty in determining DR and of the oscillatory mode. 

3. Uncertainty in determining the amplitude of the oscillatory mode. 

On the basie of the above uncertainty at high frequency in adding a 

vector representing the Fourier Transform of the oscillatory mode, it may 

be concluded that the greater the time, t^, the less the uncertainty. This 

may be observed in eq . (D-3) . To simplify that expression consider 

At = At , then 
o i 
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n 3 1 12 1 * • »p 



n = l*2 f 3 



• • • • 



i -jp^' 

[?T]q , Ai) e 
L jn o(oec; J 



-Jnw^-r 

[M][UT(t)] 1 l q 1 (T n )e 



n = 1,2 



• • • * 



The larger t becomes, the greater the percentage of the graphical analysis 



that is accomplished with the same degree of approximation in both numerator 
and denominator. If the oscillograph record were long enough, the best results 
at high frequency would be obtained by continuing summation of the output 
vectors until such time ae the end point is clearly indicated. This allows 
the same degree of approximation in obtaining both the output vector and 
the input vector throughout the entire summation. However, in addition to 
the tedium of plotting a large number of vectors, uncertainties are intro- 
duced by the oscillograph record in picking off ordinates well out in the 
time response. The noise to signal ratio becomes high, and drift of the 
baseline appears to be one of the evils of the long record. 

The most uncertain of the graphical methods for use with the highly 
oscillating type response is the method which subtracts the pure 
oscillatory mode from the entire response. In this case, the subtracted 
oscillatory mode represents a high percentage of the total response. This 
large portion of the response is represented by a single vector which is an 
exact Fourier Transform. Consequently, the output, in the extreme case where 
the response is almost entirely oscillatory, is represented exactly, and the 
input is approximated. The result is that at high frequency, the output 
reaches a certain fixed phase angle, due to the predominance of the oscillatory 



P 
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mod*, and the input continues to shift phase angle rapidly. Therefore, 
it appears that at very high frequency the oscillatory part of the output 
will eventually lead the input. The lagging phase of the remainder should 
theoretically compensate in such a manner that the vector representing the 
total performance function will have the proper phase angle. But, the 
vector representing the remainder is very small, and therefore, although 
its phase lag may be increasing in a like manner to the denominator (input), 
its contribution to the phase angle of the vector summation is very small. 
Turthermore, the additional manipulation of the time response data in 
obtaining a breakdown into components offers added uncertainty in the vector 
summation. 

Sample records are provided to show the vector summation using method 
(3) in which the oscillatory mode is removed from t = 0 to t = 00 and 
represented by an exact Fourier Transform, and method (2) in which the 
Fourier Transform of a sine wave starting at t is added to the output 
summation for t * 0 to t = t . Figures I>-9 and D-10 respectively show these 
two cases. Note that in Fig. D-B the vector representing the oscillatory 
mode is the largest of the three vectors comprising the total response 
vector. Furthermore, it has attained its maximum lagging phase angle to 
within about 5 degrees. The phase angle from the performance function at 
= 6 by this method is -49 degrees (corrected for gyro phase angle 
from Fig. D-14) . Figure I>-10 shows the results of using the same approxi- 
mation in both input and output over the greater part of the response, method 
(2) . The addition of the vector representing the sine wave summation is 
nearly negligible in its effect on the total response vector. The phase 
angle at this same frequency by this method, (2), is -87.5 degrees. 

As previously discussed, adding the eine term involves uncertainty in 
properly locating t p . This error may be reduced to a negligible size or 



eliminated by the following raeane . Select t where the output crosses the 

P 

base line and observe the magnitude of the vector resulting in comparison 

with the magnitude of the last vector in the summation from t = 0 to t = t . 

P 

If the vectors are approximately the same size, no appreciable error will 

result from an error in locating t . If the vector representing the added 

sine term is large compared to the last vector of the summation, t should 

P 

be moved to the next null point of the time response, which will decrease 

-V (DR)(T/2) 

the magnitude of the sine wave amplitude by the factor e 

It is the opinion of the authors that the added sine wave method of 

component summation is subject to the least inaccuracies, of the methods 

investigated, in determining the performance function of a response which is 

essentially oscillatory. The method affords the same degree of approximation 

in both the input and response, until such time as the response is subjected 

to inaccuracies of noise and instrumentation. At this point (t ) the small 

P 

remainder is represented by an exact trnnsform. 

Referring again to Figs. 3-3 and 3-4, the effects of various analysis 
methods in obtaining the performance function are shown. The method numbers 
are those listed in Table 2. As previously noted, the variation in PA from 
the calculated performance function at high frequencies is very high, and it 
is difficult to explain the variation in terms of the dynamics of the air- 
plane. The assumptions listed in setting up the theoretical equations of 
motion for the test vehicle, limit the phase shift to 90 degrees of lag at 
high frequencies. However, refinements which were excluded by these 
assumptions may add about 10 to 15 percent additional phase lag at high 
frequency, as shown in the theoretical performance function found in Ref. 7. 

The performance function applies to the B-26, an aircraft similar in dimensions 
and mass configuration to the B-25J . The experimental PA, at high frequencies 
are in disagreement with the calculated performance functions determined both 

for the B-25 and the B-26. However, the very low AR at the frequencies of 
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uncertain phase angles makes PA uncertainty relatively unimportant in 
determining the performance functions. 



The scattering of points in the doubtful phase angle range shows the 
effect of various methods of determining the transform of the response. It 



Compare Pig. D-8 with Fig. D~ll. This indicates that W^At has not reached 
a critical value. (Ref. 3). 

The sensitivities for this performance function were determined in the 
same manner as used in the performance function for rolling velocity response 
to input aileron responses. 

As previously noted, the calculated performance function, (eq. I>-ll), 
has a very small root in the denominator. Such a small root either positive 
o* negative would not be apparent in the time response and would be implicitly 
considered to be zero by the action of drawing in the base line. To be 
consistent for comparative purposes, the same assumption must be made with 
the calculated performance function. 

Thus, to get a valid comparison let jw - 0.0013 = Jw in the calculated 
performance function. Then by writing the performance function as a 
sensitivity multiplied by a frequency function, a sensitivity that can 
be compared with the ratio of the areas under the experimental output and 
input curves is obtained. The comparison is: 



was found by halving the At^, the phase angle was negligibly affected. 




(from run no. 4665 + 6 input) 

a 




(from run no. 4662 -4 input) 

a 



In general, the results obtained bring out several points concerning 



the application of the pulse technique as a means of determining the 



performance functions of dynamical systems. Results obtained in rolling 
velocity response for aileron inputs showed very close agreement with the 
theoretical performance function throughout the entire range of frequencies 
investigated. The comparison of time response patterns or conformations 
indicate the response pattern involving the least inaccuracies in analysis 
i s one which responds quickly, reaohing its maximum amplitude in the early 
part of the time response. Under these conditions, the end point will uaually 
be well removed from the origin, and uncertainties in added vectors represent- 
ing the Fourier Transforms of oscillatory modes, right triangles, etc. will 
have a minimum effect on the AR and PA of the performance function, until 
very high frequencies and small AR prevail. 

A potential source of error lies in establishing the true baseline for 
each response and input time record. An error in the location of this line 
produces a double error in the results in that as much false area is added 
to one side of the line as is subtracted from the other side of the line. In 
the case of records having the long tail characteristic of a lightly damped 
oscillatory mode, the envelope should be drawn in determining the base line, 
particularly if there is any drift in the center of the response record. This 
is demonstrated in Fig. 8, Ref. 3. In any case, measurements from the reference 
trace to the center-line should be made at short intervals and connected, 
rather than passing a line between two measured points at the extremities of 
the record. This corrects for warping of the record incurred in the drying 
process. If nominal care is exercised in laying out the baseline, and in 
picking off ordinates from the record, the uncertainties in reproducing 
results will be satisfactorily low. 
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FLIGHT TEST RECORDS 



FLIGHT NO. 1 

Date: 28 July 1949 Duration: 2.5 hrs 

Airplane: B-25J No. 44-30328 

Crew: 

Pilot: C. 0. Bostroa 

! 

Co-pilot: F. H. MLchaelie 

Observers: J. C. Wootton 

A. A. Hollander 

Purpose of Flight: 

Instrument shakedown 

Airplane configuration: 

Gross Weight - 26,100 

C. G. Percent MAC - 27.4 

Instrumentation installed for recording elevator position, rudder position, 
rate of roll and rate of yaw. 

Flight Procedure: 

(1) Trimmed aircraft at 10,000 ft PA and 175 mph IAS. 

(2) Observed action of oscillograph of input and response traces 
using test pulses of positive and negative aileron and rudder displacements. 

(3) Obtained records of positive and negative ailerons and 
rudder pulses . 

Results and Comments: 

Frequency regulation of the inverter was very poor. Frequency varied 
from 380 cps to 405 cps . Noise level on the oscillograph records was so 
high that the records were not satisfactory for reduction. 
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FLIGHT 50. 2 

Date: 2 August 1949 Duration: 1.5 hrs 

Airplane: B-25J 5o. 44-30328 

Crew: 

Pilot: C. Collins 

Co-pilot: J. B. Bain 

Observers* F. H. Michaelie 
J. C. Vootton 

Purpose of Flight: 

To obtain oscillograph records of airplane response to aileron and 
rudder inputs. 

Airplane configuration: 

Gross weight- 26,100 
C. G. Percent MAC - 27 .4 
Flight Procedure: 

(1) Trimmed at 10,000 ft PA and 175 mph IAS. 

(2) Observed action of oscillograph traces. for test pulses of positive 
and negative aileron pulses. 

(3) Obtain records of positive and negative aileron and rudder pulses. 
Besults and Comments: 

No records obtained due to malfunctioning of recording equipment. 
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FLIGHT HO. 3 

Date: 2 August 1949 Duration: 1.5 hrs 

Airplane: B-25J No. 44-30328 

Crew: 

Pilot: C. Collins 

Co-pilot: J. B. Bain 

Observers: F. H. Michaelis 

J . C . Wootton 

Purpose of Flight: 

To obtain oscillograph records of aircraft response to aileron and 
rudder inputs. 

Airplane configuration: 

Gross weight - 26,100 
C. G. Percent MAC - 27.4 
Procedure this Flight: 

(1) Trimmed at 10,000 ft PA and 175 rnph IAS. 

(2) Observed action of oscillograph traces for test pulses of positive 
and negative aileron pulses. 

(3) Obtain records of positive and negative aileron and rudder pulses. 
Results and Comments: 

The records obtained had an unacceptable noise level, but were somewhat 
improved over preceding flights. Input pulses, while acceptable, were too 
long for ease of data reduction. 
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FLIGHT HO. 4 

Date: 5 August 1949 Duration: 1.60 hr 

Airplane: B-25J Ho, 44-30328 

Crew: 

Pilot: Capt . Warwick 

Co-pilot: J, B. Bain 

Observers: J. C. Wootton 

F. H. Michaelis 

Purpose of Flight: 

To obtain oscillograph records of airplane response to aileron and 
rudder displacements. 

Airplane configuration: 

Grose weight - 26,100 
C.G. Percent MAC - 27 .4 

Replaced faulty inverter, and installed new Lord shock mounts on gyro can. 
Flight Procedure: 

(1) Trimmed at 10,000 ft PA and 175 mph IAS. 

(2) Observed action of oscillograph traces for test pulses of 
positive and negative aileron pulses. 

( 3 ) Obtain records of positive and negative aileron and rudder pulses. 
Results and Comments: 

Flight discontinued due to wiring difficulties. 
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FLIGHT NO. 5 

Date: 5 August 1949 Duration: 1.5 hre 

Airplane: B-25J No. 44-30328 

Crew: 

Pilot: Capt . Warwick 

Co-pilot: J. B. Bain 

Observers: J. C. Vootton 

F. H. Michaelis 

Airplane configuration: 

Gross weight - 26,100 

C. G. Percent MAC - 27.4 

Changes since last flight: 

Loose connection in 300 volt system found and repaired. 

Flight Procedure: 

(1) Trimmed at 10,000 ft PA and 175 rnph IAS. 

(2) Observed action of oscillograph traces for test pulses of 
positive and negative aileron pulses. 

(3) Obtain records of positive and negative aileron and rudder pulses. 
Results and Comments: 

The records obtained were considered to be useable. The noise was of 
insignificant amplitude, and the input pulses were of the correct size, 
however, the roll gyro picked up a vibration, apparently from its mount, 
signal of about 10 cps which would make data reduction somewhat difficult. 
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FLIGHT NO. 6 

Date: 9 August 1949 Duration: 1.5 hrs 

Airplane: B-25J No. 44-30328 

Crew: 

Pilot : C . Collins 

Co-pilot: J. B. Bain 

Observers: F. H. Michaelis 

F. Smith 



Airplane Configuration: 

Grose weight - 26,100 

C. G. Percent MAC - 27.4 
Changes since last flight: 

Intailed stop to remove shock mounts while airborne and wedged gyro 
container rigidly to airplane to eliminate vibration of the mount from 
being recorded. 

Procedure this flight: 

(1) Trimmed at 10,000 ft PA and 175 mph IAS. 

(2) Observed action of oscillograph traces for test pulses of positive 
and negative aileron pulses. 

(3) Obtain records of positive and negative aileron and rudder pulses. 
Results and Comments: 

The records obtained on this flight were very good. The attenuator 
settings used were: 

Function Attenuator 



"z 3 

W x 4 

6 5 

a 

6 r 6 

These attenuations allowed the traces to remain on the pa£er and gave large 

enough traces to be readable. Eight records were taken and static 
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calibrations with the master pickoff were made in flight. 
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FIGURE D-l DETERMINATION OF PERFORMANCE FUNCTION, [PF] (A)(&<lVV|)l AT W f -1.25 RAD/SEC BY METHOD 

OF ANALYSIS (1) OF TABLE 2 
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FIGURE D-2 DETERMINATION OF PERFORMANCE FUNCTION, [PF ] (A)(4oV *), AT W f -1.25 RAD/SEC BY METHOD 

OF ANALYSIS (2) OF TABLE 2. 
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FIGURE D-4 DEMONSTRATION OF USE OF METHOD OF ANALYSIS (2) OF TABLE 
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FIGURE 07 DETERMINATION OF PERFORMANCE FUNCTiON. CPFl,^ , AT THE RESONANT FREQUENCY OY METHOD NO. 
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FIGURE D-d. DETERMINATION OF PERFORMANCE FUNCTION, CPF ] (A ) (8aWl ) , AT W f = 5 RAD/ SEC BY METHOD 

OF ANALYSIS (2) OF TABLE 2. 
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FIGURE D -10. DETERMINATION OF PERFORMANCE FUNCTION, [PF] {A)c8 Wz) AT W f « 6 RADS/SEC BY 

METHOD OF ANALYSIS C2) OF TABLE 2. 
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FIGURE D ~ 1 1 . DETERMINATION OF THE PERFORMANCE FUNCTION, CPF] (a)(5qW;) , FORW-S RADS/SEC 
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MASSACHUSETTS INSTITUTE OF TECHNOLOGY 
AERONAUTICAL ENGINEERING DEPARTMENT 
INSTRUMENT SECTION 



(ARSR) 
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APPEND II £ 



GLOSSARY 



The system of notation adopted for use in this thesis has been 
selected because of its peculiar usefulness in representing operating 
system performance. In particular, it has the following specific 
advantages : 



1. It is easily learned. 

2. It is adopted to a wide range of situations. 

3. It is built up almost exclusively of characters found on the 
keyboard of a standard American typewriter. 

4. Any one of the compound symbols of the system is readily inter- 
preted without recourse to an extensive glossary. 

A representative list of the primary symbols and those used as sub- 
scripts is given in the following table. 

Primary Symbole 



(AH) 



( )(< »i n Vt> 



(CT) 

(DR) 

(fr) 

(FT) 

s 



Angle 

Input - output amplitude ratio for a 
given operating component = 

q (in)a 

Linear viscous damping coefficient 

Characteristic time 

Damping ratio = - 

2 

Frequency ratio = 

n 

Fourier transform 

Force of gravity/unit of mass 

Moment of inertia of airplane about 
I axis 





Product of inertia 
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1 zz 


Moment of Inertia of airplane about 
Z axle 


L 


Bolling moment 


L 

P 


Bolling moment due to rolling velocity 


L 

r 


Bolling moment due to yawing velocity 


L * a 


Bolling moment due to aileron dieplacement 


L 

T 


Bolling moment due to side slip velocity 


(LT)[q(t)] 


Direct Laplace transform of q(t) = Q(s) 


m 


Mass 


N 


Tawing moment 


N 

P 


Tawing moment due to rolling velocity 


N 

r 


Tawing moment due to yawing velocity 


N 

V 


Tawing moment due to side 9lip velocity 


N 6 8 


Tawing moment due to aileron deflection 


"<>r 


Tawing moment due to rudder deflection 


^ ><«in Vt> 


Input - output performance operator for 
a given operating component 


(PA> ' >S n Vt> 


Input - output dynamic phase angle for a 
given operating component 


(P?)( >Sn Vt> 


Input - output performance function for a 
given operating component, where 

j[ V *(pa) ( i b ,] 

q ( in) q (in)a 6 

* (PA) (out)] 

q (out) = q (out)a * 


q (in) 


Input quantity 


q (out) 

q(t) 


Output quantity 

A quantity which is a function of the 
real variable, time 


BT 


Bight triangle 


8 


Complex variable used in Laplace transform 
theory 


S < >S B 1.at> 


Static sensitivity for a given operating 
component 
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*d 


Time delay 


At 


finite increment of time 


U 


Airplane forward velocity 


(UT) 


Unit triangle 


(URT) 


Unit right triangle 


V 


Acceleration along T axis 


w ( ) 


Airplane angular velocity about axis ( ) 


w f 


Angular forcing frequency 


W n 


Natural frequency 


T 


Side force 


T 

V 


Side force due to side slip velocity 



A dot over a variable represents its derivative with respect to time 

(e .g. ( ) = “ ^ ) . 

Two dote over a variable represents a second derivative with respect to 

time (e.g. ( ) « - — ). 

dt^ 

The following Greek letter symbols are retained because their meanings have 
become thoroughly identified with these symbols in the aeronautical engineering 



field . 




6 

a 


Total aileron deflection (sum of left and 
right aileron deflection) (positive deflection 
of right aileron is downward) 


6 

r 


Rudder deflection (left rudder deflection is 
posit ive) 


* 


Angle of bank 




Angle of yaw 


u> 


Angular forcing frequency 
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Modifying symbols (subscripts) 



a 


Aileron 


(*) 


Amplitude 


A 


Airplane 


(dir) 


Divergent 


(in) 


Input 


(ind eye) 


Indicating system 


M 


Measuring 


0 


Output 


(out) 


Output 


(08C ) 


Oscillating 


P 


Bate of change of bank angle ~ 


(rem) 


Remainder 


r 


Rate of change of heading = V 

L 


▼ 


Side slip velocity 
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